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INVESTIGATION  OF  A  COMPLEX  TECHNIQUE  OF  SMOKE 
PARTICLE  DEPOSITION  ON  SCAVENGERS 


1.  INTRODUCTION 

The  research  program  "Investigation  of  a  Complex  Technique 
for  Smoke  Particle  Deposition  on  Scavengers  (SPADES)", 
DAAK-11-83-K-0007  concentrated  on  several  promising  results  of 
the  previous  investigations  under  DAAK-11-81-C-0075  and 
DAAG-29-79-C-007 3 .  It  stressed  the  deposition  of  smoke 
particulates  on  nonspherical  bodies  which  had  shapes  of  planar 
scavengers  (disks,  plates,  and  grids),  falling  steadily  or 
performing  an  oscillatory  motion.  Special  attention  was  oaid  to 
the  aerosol  collection  efficiency  of  scavengers  made  of  fine 
electrete  fibers  organized  in  a  regular  mesh  or  forming  a  thick 
plate  of  randomly  oriented  fibers  (Dust  Magnet,  Poly  Mag  80  and 
90,  product  of  Brudhow  Industries,  Collierville,  TN).  These 
experiments  were  aimed  at  confronting  the  results  of  the 
numerical  modeling  of  smoke  particle  deposition  on  uncharged  and 
charged  collectors  with  the  measurement  in  a  laboratory  wind 
tunnel . 

Most  of  the  laboratory  experiments  have  been  performed  with 
smoke  particles  generated  by  the  reaction  of  titanium  chloride 
with,  water  vapor.  This  reaction  with  the  subseguent  particle 
shattering  in  an  electric  furnace  yielded  a  narrow  size 
distribution  of  smoxe  particles  with  the  modus  lyinq  within  the 
size  range  of  "Greenfield  gap"  (0.2  pm  <  d  <  0.8  urn).  The 
technique  of  aerosol  generation  has  been  described  in  the  Final 
Report  DAAX-1 1-8 l-C-0075 .  Our  main  task  under  this  grant  was  to 
increase  the  reproducibility  and  duration  of  smoke  particle 
generation  for  the  investigation  of  particle  deposition  on  models 
in  a  wind  tunnel,  or  for  the  light  extinction  measurements  in  a 
smoke  chamber. 

One  of  the  main  tasks  in  the  proqram  DAAK-11-8 3-K-0007  was 
the  investigation  of  the  behavior  of  a  cloud  of  falling 
scavengers  in  which  the  collectors  do  not  interact 
hydrodynamical ly ,  however,  as  a  whole  they  induce  a  downdraft. 
That  on  one  side  might  entrain  the  smoke  particles  and  transport 
them  downwards,  however,  on  the  other  side  the  decreasina 
relative  velocity  of  a  smoke  particle  and  collector  leads  to  the 
lowering  of  the  collection  efficiency  of  the  scavenger. 

Finally,  an  assessment  of  the  total  mass  (number  of 
scavengers  of  a  specific  type  for  improving  the  visibility  in  a 
smoke  cloud  was  made.  This  assessment  accompanied  by  a  suggested 
scavenger  dispersion  over  a  smoke  zone  gives  an  idea  about  the 
practical  applicability  of  the  SPADES  research  proararn  results. 
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2.  SCAVENGER  COLLECTION  EFFICIENCY-THEORETICAL  STUDY 


The  premise  of  this  program  was  that  a  planar  collector  ha9 
a  higher  scavenging  efficiency  than  a  sphere  (drop)  of  the  same 
cross  sectional  area  (Martin  and  Podzimek,  1982)  and  the  same 
volume  of  displaced  fluid.  Because  of  the  given  size  of  smoke 
particles  (0.2  <  d  <  1.0  um) ,  the  original  model  used  by  the 
investigators  at  UCLA  (Pitter,  1977;  Martin  et  al.,  1980a,  b,  c) 
had  to  be  modified  and  extended  into  the  large  Reynolds  numbers 
and  checked  the  suitability  and  additivity  of  terms  describing 
the  particle  deposition  by  Brownian  diffusion,  phoretical  and 
electrostatic  forces. 


Assuming  that  the  mass  of  smoke  particle,  m  ,  is  much 
smaller  than  the  mass  of  the  collector  (e.q.,  disk,  grid)  one  can 
deduce  the  aerosol  particle  trajectory  from  the  equation  for  the 
particle  velocity,  v  , 


d  v 

-k*  =  m  5 

dt  p  ’ 


6  it  vi  r 

a  r%  -¥  -*■ 

> . -  ,  (V  -  U  )  +  F_.  +  F_  .  ,  +  F_ 

(1+aKn)  pa  Th  Dlf  B 


(1) 


This  equation  does  not  include  the  term  corresponding  to  the 
effect  of  Brownian  diffusion,  which  is  applicable  to  particle 
radii,  r  <  0.1  um.  In  a  separate  numerical  model  it  was 
includedpearlier  together  with  the  effect  of  electrostatic  charge 
by  Martin  et  al.  (1980c).  However,  the  terms  for  Brownian 
diffusion,  thermo-  and  dif fusiophoretical  forces  can  be  neglected 
in  comparison  to  the  effect  of  inertial  and  electric  forces  for 
the  TiCl.  smoke  particles  of  the  size  distribution  curve  analysed 
and  published  earlier  (Podzimek,  1983).  For  the  same  reason  the 
slip  flow  correction  term  a  Kn  witha  ■  1.26  +  0.44  x  exp  (-1.10 
kn-1)  and  Kn  *  can  be  neglected. 

rP 

From  the  known  particle  trajectory  around  the  disk  (thin 
oblate  spheroid  with  axes  ratio  of  0.05)  the  disk  collision 
efficiency 


E 


*(W 


(2) 


was  deduced.  Yc  is  the  largest  offset  of  the  particle  from  the 
disks  main  axiscthe  particle  can  have  and  still  collide  with  the 
crystal.  No  attempt  was  made  to  speculate  about  the  potential 
"annular  effect"  decreasing  the  number  of  deposited  particles 
around  the  stagnation  point  on  the  obverse  side  of  the  disk 
(Pitter  and  Pruppacher,  1974). 


A  composite  diagram  for  different  Re  (up  to  Re  =*  80), 
uncharged  particulates  is  presented  in  Fig.  1.  It  assumes  the 
validity  of  Eg.  (1)  in  which  phoretical  terms  <FTh/  Fn'f  '  and 
Brownian  diffusion  are  small  in  comparison  to  theinertill 
deposition  for  particle  radii  greater  than  0.1  um.  The  formulas 
for  creeping  flow  conditions  and  for  the  forces  acting  on  an 
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aerosol  particle  were  described  elsewhere  (Martin  et  al,  1980a, 
l°80b,  1980c;  Martin  and  Podzimek,  1982;  Podzimek,  1983).  The 
most  significant  result  one  can  deduce  from  Pig.  1  based  on  the 
theoretical  calculations  for  a  disk  is  that  the  collection 
efficiency  for  smoke  aerosol  particles  (with  radii  between  0.1  pm 
and  2.0  /am)  will  not  surpass  8%.  The  inertial  deposition  will  be 
the  main  mechanism  for  smoke  particle  scavenging  with  a 
considerable  influence  of  thermophoretical  forces  (if  they  exist) 
in  the  particle  size  range  0.3  <  r  <  0.8  ym.  The  only 
explanation  of  the  peculiar  behavior  of  particles  at  Re  *>  0.1; 

0.5  and  1.0  (which  shows  E  <  10”  for  all  particle  sizes 
considered  and  therefore  is  not  drawn  in  Pig.  1)  is  the  basic 
change  of  the  flow  pattern  around  the  oblate  spheroid  at  low  Re 
and  the  evolution  of  vortices  on  the  reverse  side  at  Re  >  5.0. 

At  Re  <  0.5  particle  laden  air  is  longer  time  in  contact  with  the 
collectors  surface. 


Much  attention  has  been  paid  to  the  analysis  of  different 
kind  of  electric  charging  of  the  collector  and  smoke  particle  as 
well  and  to  its  effect  on  collection  efficiency  of  a  planar 
collector  (Podzimek  and  Martin,  1984).  In  spite  of  the  fact  that 
both  field  and  diffusion  charging  might  affect  the  particle 
deposition  a  simplified  model  was  applied  to  the  force  actinq  on 
aerosol  particle  in  the  proximity  of  the  scavenger 


Fl,2  *  "19 


♦  - 


6narl 
(1+2  Kn) 


(Vrv2)  +  7T  +  K27pv  +  01E2  (3) 


and  on  the 
particle 

fu 


collector  due  to  the  image  charge  induced  by  the  smoke 


m 


2g*  • 


6-n-n  r_ 

Cl  4 


<vv 


VT  + 


K27ov 


+  0- 


2  &1 
corrected 


4  ) 


(l+2aKn) 

The  terms  in  each  equation  represent  the  buoyancy 
gravitational  force,  corrected  drag  in  accordance  with  Stokes 
resistance  law  (assumption  is  made  that  the  effect  of  planar  and 
spherical  collector  are  similar).  Further  it  was  assumed  that 
the  scavenger  position  and  orientation  of  the  collector's  axes 
remains  undisturbed  by  either  hydrodynamic  interaction  of 
particles  or  by  electrical  forces  (Martin  et  al,  1981).  The 
Coulomb  forces  are  the  dominant  ones  and  are  radiatinq  from  the 
mass  center  of  each  body  (image  forces  are  not  considered).  The 
justification  of  creeping  flow  equations  and  the  numerical 
technique  for  a  thin  oblate  spheroid  has  been  discussed  by  Martin 
et  ai.  (1980b).  Under  this  program  additional  flow  fields  for  Re 
=  80  were  developed.  There  were  considerable  difficulties  with  a 
similar  model  for  Re  =  100  due  to  the  flow  instability.  In  case 
that  there  is  no  evaporation  or  vapor  condensation  on  smoke 
particles,  the  d i f f us iophoret ic  and  thermophoret ic  terms  were  set 
to  zero.  The  collection  efficiency  curve  for  Re  =  80  is  compared 
to  the  old  curves  (Martin  and  Podzimek,  1982)  for  low  Re  and 
plotted  in  Fig.  2.  In  conclusion,  the  electric  charge  effect  can 
dramatically  change  the  particle  deposition  on  a  charged 


10 


-If  / 


Re  •  30^-." 


n  >'  *  * 
0.1  -/  • 

/o 


/ 


**  /  , 

>  i 
r  i  i 
/  /  i 

'  /  i 

/  ; ! 

*  /  / 


R«*50  ✓  ' 


0.03 


.<*  6  1, _ . 

r  *'•*  "*■* 

'*  *  9 

I  •  • 


/  / 
/  / 


1  10,' 
eT”*- _ - 

V  &  -  -  0.1 

/  •<  *v 


/r*t 

IM\\ 

a0,JJ;  »\ 

h  >  \ 

!/  .  V 


\  \s.O 
'  \ 

2.0V.  ' 


a  2  kv 
*5  kV 

o3  kv  for  3sec. 
•  Q  kV 

O' rim  deposit 
O i  kV  w/  Au  coat 
*5  kV  M  “  ' 

a  no  el.  charge 


5  fp  (/im) 


Fig.  2  Collision  efficiency  of  charged  and  uncharged  oblate 
spheroid  (disk)  as  a  function  of  Re,  particle  radius, 
r  ,  and  electric  charge. 
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collector.  The  surface  charge  assumed  in  our  calculation  is  high 
(0  *  2.0  esu  cm”2)  indeed  (2.0  esu  cm”2)  and  it  can  be  hardly 

maintained  in  a  dense  aerosol  cloud.  At  the  beginning  of  the 
particle  deposition  it  might  enhance  the  particle  deposition  of 
one  order  of  magnitude.  The  other  interesting  feature  of  the 
assumed  positive  interaction  of  smoke  particles  with  a  charged 
collector  is  the  inefficiency  of  such  a  process  for  particulates 
with  radii  larger  than  5.0  urn  if  compared  to  the  collision  of  not 
electrically  charged  particles.  This  is  also  apparent  in  Fig.  3 
where  is  plotted  the  calculated  collision  efficiency  of  a  thin 
oblate  spheroid  as  a  function  of  Stokes  number.  K. 

Considerable  effort  has  been  made  to  investigate  the 
collision  efficiency  of  cylinders  (fibers).  The  collision 
efficiencies  calculated  from  models  presented  by  different 
authors  are  plotted  in  Fig.  4.  There  is  still  a  considerable 
uncertainty  of  different  models  in  the  domain  of  low  Re,  however, 
it  can  be  shown  that  our  data  (curve  4  in  Fig.  1)  for  potential 
flow  are  not  much  different  from  the  calculation  by  Davies  and 
Peetz  or  from  the  stepwise  calculation  by  Harrop  (using  Happel 
flow  field)  -  curve  1  in  Fig.  4.  The  calculation  of  the  flow 
field  and  particle  trajectories  are  presented  in  Appendix  A. 

There  is  still  a  gap  in  our  knowledge  of  the  ultrafine  particle 
deposition  in  a  viscous  flow  around  a  very  fine  fiber,  what  would 
mean  to  extend  our  curves  in  Fig.  4  into  the  domain  of  Stokes 
numbers  K£  0.01.  For  instance,  the  Kuwabara  model— used  by 
several  authors  and  supported  roughly  by  experiments  (Kirsch  and 
Fuchs,  1967)  for  the  determination  of  particle  trajectories 
around  a  fiber — is  too  sensitive  to  the  starting  conditions.  The 
same  holds  for  the  model  used  by  Schlamp  (1977)  for  calculation 
of  cloud  droplet  deposition  on  columnar  type  crystals,  which  the 
author  approximated  by  infinite  cylinder  (Appendix  B). 

The  main  results  of  the  investigation  (Liu,  1985)  for  ideal 
flow  without  electrostatic  force  and  aerosol  particles  negligibly 
small  compared  with  the  cylinder  can  be  summarized  as  follows: 
There  is  a  reasonably  good  agreement  between  our  data  and  the 
collision  efficiency  calculated  by  Albrecht  (1931),  Sell  (1931) 
and  Glauert  (1940),  Langmuir  and  Blodgett  (1946)  and  Davies  and 
Peetz  (1955)  for  the  cylinder  in  the  transition  flow  region  (Fig. 
5).  The  slight  disagreement  of  the  curves  can  be  explained  bv 
the  higher  sensitivity  of  the  numerical  calculation  with  respect 
to  the  starting  point  ahead  of  the  cylinder.  Glauert  and  Sell 
and  Albrecht  started  at  three  cylinder  diameters  where  they 
assumed  aerosol  particles  having  the  same  velocity  like  an 
airflow  at  infinity.  Langmuir  and  Blodgett  started  at  four  and 
Davies  and  Peetz  at  five  diameters.  In  our  case  we  selected  six 
diameters  and  hoped  to  obtain  sufficiently  accurate  particle 
trajectory  along  the  obverse  side  of  a  cylinder.  Another  reason 
for  the  disagreement  of  curves  in  Fig.  5  might  be  the  different 
particle  density  assumed  by  the  investigators.  In  Fig.  6  are 
plotted  the  particle  trajectories  about  a  cylinder  positively 
charged,  uncharged  and  negatively  charged.  This  causes  a  5.0  ^m 
(radius)  particle  to  be  repulsed  or  attracted  close  to  the 
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Fig . 


3  Calculated  collision  efficiency  of 
a  thin  oblate  spheroid  as  a  function 
of  Re  and  Stokes  number,  K. 


4  Calculated  collision  efficiency,  E, 
of  an  infinite  cylinder  as  a  function 
of  Stokes  number  of  particles,  K,  by 
different  authors. 
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Fig.  5  Comparison  of  theoretical  collision  efficiencies 
for  a  circular  cylinder  and  aerosol  particles  as 
a  function  of  Stokes  number. 


surface  of  a  metallic  cylinder  in  addition  to  the  action  of 
inertial  forces.  The  particle  density,  p  »  1.0  q  cm”  in  this 
case,  plays  an  important  role  in  modifying  the  trajectory  and  the 
number  of  deposited  particles,  as  it  is  demonstrated  in  Fin.  7, 
where  for  comparison,  the  particle  trajectories  are  plotted  for 
p  ■  2.0  g  cm  .  The  effect  of  positively  and  negatively  charqed 
collector  on  the  aerosol  particle  (of  the  same  size  mass  and 
charge)  collision  is  different  (in  Pig.  8  are  plotted  the 
collision  efficiencies  for  the  example  case  described  in  Fig.  6). 
The  comparison  of  Figs.  6  and  7  leads  also  to  the  conclusion  that 
for  enough  strong  attractive  Coulomb  force  the  collision 
efficiency  will  be  qetting  larger  for  smaller  particle  density 
( p  ■  1  g  cm-3).  For  repulsive  Coulomb  force  the  situation 
reverses  if  the  electric  charge  force  overrides  that  of  the 
inertial  force. 

For  the  sake  of  estimating  the  potential  effect  of  the 
h'/drcdynamic  interaction  of  two  cylinders  in  the  transition  flow 
region,  we  studied  the  particle  deposition  on  two  cylinders  close 
together.  The  main  results  of  these  theoretical  stud les— which 
are  described  in  more  detail  in  Appendix  C-- -can  be  summarized  as 
follows:  Collision  efficiencies  for  an  example  case  (D  »  0.1 
cm,  U  •  100  cm  a”1,  P  ■  1.0  g  cm”3)  of  particle  deposition  in 
potential  flow  are  increasing  with  diminishing  distance  between 
cylinder  (fiber)  centers,  S,  until  approximately^—  ■  2.0,  and 
than  decrease  with  smaller  value  of£/Dc  (Fig.  9).c  At  a  value  S/D 
*  1.2  the  two  parallel  fibers  act  approximately  like  a  single 
body  (Shyu,  1986). 

An  attempt  was  made  to  substitute  a  rather  unrealistic 
application  of  a  potential  flow  to  the  not  very  close  packei 
parallel  fibers  by  a  model  simulating  an  airflow  through  a  many 
fiber  array  forming  a  rectanqular  net  (channel  model)  or  a 
staggered  model  (Brown,  1984).  These  models  are  substituting  The 
stream  functions  suggested  by  Kuwabara  (1959),  Happel  (1959)  cr 
Spielman  and  Goren  (1969)  of  the  type 

*  -  IT  ll'l+2  llnl>  ,tne  '  <c' 


where  (  is  a  specific  hydrodynamic  factor,  by  many  -  fiber  term 

*  •  *  V  '  k^0  v  ,in  t1  e0B  *ir  ' ' 


The  many-fiber  stream  function  reflects  (in  the  form  of  double 
Fourier  series  in  x  and  y)  the  spacial  periodicity  ad  might  be 
suitable  for  the  investigation  of  smoke  particle  deposition  on 
filter  type  electrete  fiber  scavengers.  The  meaning  of  the 
geometrical  parameters  i  and  e  is  apparent  from  the  description 
in  Appendix  D.  In  conclusion,  it  seems  suitable  to  our  task  of 
Investigating  the  deposition  of  smoke  particulates  on  fiber  type 
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Fig.  6c  Particle  trajectory  with  repulsive  force,  p  *  1.0  g/cm 

Dc  *  0.05  cm,  Dp  «*  10  pm,  N{^_  «  55,  AR  *  1?6,  QA  *  7.5  x 
lO-2  e.s.u.  (Qav  -  20),  Qb  =  5.0  x  10“6  e.s.u.  (Qq *  -  20) 


Fig.  7a  Particle  trajectory  without  electrostatic  force, 

2.0  g/cm^,  D  ■  0.5  cm,  =  10  pm,  N_  =  55,  AR 
c  p 


Fig.  7b  Particle  trajectory  with  attractive  force,  p  =  2.0  g/cm 
Dc  =  0.05  cm,  D  =»  10  um,  NRe  *  55,  AR  *  1/6?  QA  =  7.5  x 
10"2  e.s.u.  (Qav  »  20),  Qb  »  -5.0  x  10  6  e.s.u.  (QA*  *  - 


Particle  trajectory  with  repulsive  force,  Op  -  2.0  g/cm 
Dc  =  0.05  cm,  Dp  =  10  um,  Nfce  =  55,  AR  =  1/6,  Qa  -  7.5 
10"2  e.s.u.  (Qa*  *  20),  Q_  =  5.0  x  10“6  e.s.u.  (Q  '  *  2 


Fig.  8  Numerical  collision  efficiency  E  as  a  function 
of  particle  Stokes  numbers  K  and  Dc  =>  0.05  cm, 
NRe  »  55,  o  ■  1  g/cm3  and  AR  ®  1/6. 
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scavengers,  to  develop  further  Brown’s  model  because  the  model  is 
relatively  simple  and  it  gives  a  reasonably  good  description  of 
the  flowfield.  The  last  statement  is  based  on  Brown's  analysis 
and  on  our  preliminary  experience  which  showed  also  a  limited 
applicability  of  the  model  for  thick  filter  material  (pressure 
drop  and  changing  Re  across  the  filter).  The  inaccuracy  of  the 
model  very  close  to  the  fiber  surface  does  not  affect  much  the 
deposition  of  particles  the  si2e  of  which  excludes  a  strong 
influence  of  diffusional  deposition.  A  more  complicated 
expression  for  stream  function  will  be  required  for  a  staggered 
model  of  a  filter  (Brown,  1984), 
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which  with  corresponding  boundary  values  will  enable  to  calculate 
the  flowfield  in  the  array  of  fibers  described  in  Appendix  D. 

The  rate  of  smoke  particles  deposited  on  each  row  plan 
perpendicular  to  airflow  direction  can  be  calculated  according  to 
the  formula 


dni  (rpi,t) 


-2R  L  Np^VRjni  (  rpi  RirlPp* 


(8) 


which,  after  integration,  yields 


nj  (rpi't)  *  nj  (rpi'0)  exp 
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Vl  (rpi'C)  *  (  i“n  j  (  Tpi  ,  t )  !  exp  (-  Aij  +  1t) 

and  for  a  constant  A.  .  =  A..  ^  =  ^ii+2  S  ^  and  the  same  time 
interval  t.,  the  number  of^moke  particulates  remaining  in  the 
airflow  after  passing  e.g.,  four  planes  of  fibers  will  be 

nrem  *rpi,4t*  *  n ( rp£ , 0 ) ( 1-exp ( -A , t )  +  exp(-2A  .t)  -exp  {  —  3 A  .t) 

+  exp  ( —  4 A . t ) ]  -  3  exp  ( —  A . t )  +  2  exp  (-2A.t)  -  exp  (-3A.t)  .  (9) 


This  model,  which  assumes  that  there  is  not  a  considerable 
pressure  drop  on  individual  fiber  planes  seems  to  be  very 
appropriate  to  any  type  of  filter  or  qrid  type  scavenger  where 
the  assumption  of  parallel  fibers  is  acceptable.  The  problem  of 
grid  type  scavengers  with  node  points  or  the  effect  of  randomly 
oriented  fibers  requires  further  study. 
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3.  EXPERIMENTAL  INVESTIGATION  OF  THE  CLEARING  OF  A  SMOKE 
CLOUD  BY  "SPADES"  TECHNIQUE 

Smoke  particle  deposition  on  scavengers  (SPADES)  was 
investigated  mainly  in  the  laboratory  wind  tunnel  and  in  the 
cylindrical  smoke  chamber  (1.6  m  high  and  0.8  m  in  diameter). 

The  aim  of  these  experiments  was  the  following:  To  complete  the 
studies  made  under  the  U.S.  Army  Research  Office  Grants  DAAG 
29-79-C-0075  and  DAAK  1 1-8 l-C-0075 ,  and  to  concentrate  mainly  on 
grid  type  scavengers  and  planar  scavengers  made  of  electrete  tvoe 
fibers.  These  fibers  were  much  more  effective  than  solid  planar 
scavengers  (e.g.,  disk),  the  collection  efficiency  of  which  was 
about  one  or  two  percent  for  TiCl.  smoke  particles  with  modus 
diameter  around  0.4  um.  An  oscillatory  motion  did  enhance  only 
slightly  the  collection  efficiency  of  thin  disks  (Podzimek, 

1983). 


Electrete  type  fibers  showed  high  collection  efficiency  (7% 
for  particle  diameters  smaller  than  0.5  ^m  at  an  airflow  velocity 
of  0.25  m  s”  .  This  indicates  the  important  role  played  by  the 
diffusion  and  electrostatic  charge.  The  additional  charging  of 
electrete  fibers  will  not  enhance  dramatically  the  deposition  of 
smoke  particles  with  diameters  smaller  than  l.OjUm  (Fig.  10). 

For  the  same  electrete  fiber  diameter  (0.004  cm)  a  higher  flow 
velocity  (1  m  s“A)  stressed  a  little  more  the  effect  of  the 
additional  electrostatic  charge,  however,  the  total  effect  was 
smaller  than  in  the  previous  case— in  maximum  6.5%  (5%  for 
uncharged  fiber)  for  particle  diameters  smaller  than  0.5  (Fici. 
11).  The  other  very  important  effect  is  the  chanqinq  of  size 
distribution  curve  of  smoke  particles  which,  after  the  first 
dropping  of  electrete  scavengers,  will  remain  in  the  air. 
Scavenger  fibers  analyzed  under  a  scanning  electron  microscope 
reveal  that  size  spectrum  of  particles  found  on  the  fiber  is 
shifted  towards  larger  sizes,  and  in  consequence  of  that,  in  the 
air  will  remain  smaller  particles  (Fiq.  12).  This  result  could 
be  anticipated,  however,  what  is  the  most  important,  that  this 
shifting  in  smoke  particle  size  distribution  would  reauire  the 
use  of  different  scavengers  (e.q.,  the  use  of  fast  falling  coarse 
fiber  scavengers  first  and  then  dropping  of  finer  electrete 
fibers  of  small  settling  velocity  enhancinq  more  deposition  of 
fine  particles  by  diffusion).  The  other  possibility  of  using  the 
most  effective  scavengers  is  to  combine  planar  coarse  fiber 
scavengers  with  fine  fiber  net  stretched  over  the  coarse  fibers 
ir  order  to  remove  the  larqest  portion  of  the  aerosol  size 
spectrum  (on  this  point  a  Patent  Application  will  be  submitted 
soon,  through  the  U.S.  Army  authorities). 

The  very  important  part  of  experiments  was  the  assessment  of 
the  reproducibility  of  the  measured  collection  efficiency  of  mesh 
tyDe  scavengers  in  the  laboratory  wind  tunnel  and  in  the  0.8  mA 
smoke  chamber.  The  first  task  was  closely  related  to  the 


Standardization  of  the  TiCl^  smoke  particle  generation--which  is 
now  controlled  during  each  experiment  in  the  wind  tunnel  within 
an  error  of  approximately  +  20%  of  concentration  in  the  main 
particle  size  groups.  Thi¥  task  prompted  new  series  of 
measurements  which  led  to  the  special  one-year  research  program 
supported  by  the  University  of  Missouri  the  goal  of  which  is  the 
use  of  neutron  activation  analysis  for  determining  the  number  and 
total  nass  of  deposited  Ti-containing  particles  on  electrete 
fibers.  The  fitV:  pilot  experiments  done  in  1986  partly  within 
the  DAAK-11-8  3-K-0007  program  were  very  promising.  In  general  is 
expected  that  the  neutron  activation  performed  in  the  reactor  of 
UMR  will  enable  to  detect  small  guantities  oi  Ti  containing 
particles  down  to  0.10  pc,  with  a  mean  error  of  around  +  10%. 

This  guantity  combined  with  the  use  of  laser  cavity  aerosol 
spectrometer  for  particle  size  distribution  measurement 
represents  a  substantial  improvement  of  the  techniaue  for  the 
determination  of  the  scavenger  collection  efficiency.  The  final 
report  on  the  new  measuring  technigue  of  smoke  particle 
deposition  on  the  fibers  (UMR  grant)  will  be  sent  to  the  U.S. 

Army  Research  Office  during  the  year  1987. 

More  difficult  was  the  assessment  of  the  accuracy  and 
reproducibility  of  measuring  the  visibility  in  a  smoke  chamber 
after  a  specific  number  of  scavengers  was  dropped  into  a  dense 
TiCl.  smoke  cloud.  Besides  the  smoke  particle  generation 
technigue  (described  in  the  DAAK  11-81-C-0075  and  DA AG 
29-79-C-007 3  reports)  the  following  factors  affect  the 
measurement  of  visibility  (instrumentation  has  been  described 
also  in  the  previous  reports)  in  the  smoke  chamber:  The  type  of 
scavenger,  the  art  of  scavenger  dropping  (freguency  of  individual 
droppings,  dropping  of  scavenger  groups  in  specific  time 
intervals),  width  of  the  light  beam  for  visibility  measurement 
(probability  of  a  scavenger  crossing  the  light  beam),  size  of  the 
smoke  chamber  (air  or  smoke  circulation  induced  by  the  falling 
scavengers)  and  the  model  used  for  analysis  and  interpretation  of 
the  measured  visual  range. 

Most  of  the  experiments  have  been  done  by  the  electrete 
denoted  as  "Dust  Magnet"  Poly-Mag  80  and  90,  a  product  of  Bruhow 
Industries,  Collierville,  TN.  Several  studies  have  been  done 
with  paper  disks  and  sguares  described  elsewhere  (Podzimek, 

1981).  The  coarse  electrete  fibers  200-300  ym  in  diameter  were 
arranged  in  a  grid  type  array  2.5  cm  x  2.5  cm  (or  1.25  cm  x  1.25 
cm)  in  size  in  which  fibers  of  a  mean  diameter  of  250  um  are 
alternating  with  gaps  of  approximately  500  to  750  um.  The  mean 
scavenger  mass  was  0.105  g  and  the  fall  velocity  of  176  cm  s”  . 
The  filter  type  scavenger  had  dimensions  2.5  cm  x  2.5  cm,  x  1.3 
cm  (or  1.25  cm  x  1.25  cm  x  0.7  cm).  The  dense  packed  fine  fibers 
(approximately  40  um  in  diameter)  were  randomly  oriented  and  the 
scavenger  had  a  settling  velocity  similar  to  the  coarse  one.  The 
main  aim  of  these  experiments  was  to  estimate  the  total  clearing 
effect  in  a  dense  smoke  cloud  after  a  specific  number  of 
scavengers  was  released  and  to  relate  these  observations  to  the 
measurement  in  a  wind  tunnel.  Another  important  goal  of  these 
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Fig.  12  Size  spectrum  of  TiCl.  particles:  a  -  deposited  on  the 
fiber,  b  -  in  the  airflow. 
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Change  in  light  extinction  in  a  smoke  cloud  (Aext/min  5 
mV/min)  was  caused  by  dropping  a  specific  number  of 
scavengers  in  5  minute  intervals.  Triangles  correspond 
to  fine  fiber  "Dust  Magnet"  scavengers,  dots  to  the 
coarse  fibers. 


investigations  was  to  find  the  optimal  groups  of  scavengers  with 
the  highest  scavenging  efficiency  and  to  find  the  most  suitable 
time  intervals  between  the  dropping  of  scavengers.  Typical  size 
distribution  curve  (cumulative  size  spectrum  curve)  is  presented 
in  Fig.  12.  There  is  clearly  a  broken  line  concave 
upwards— which  supports  strongly  the  observations  made  in  a  wind 
tunnel  earlier. 

All  data  from  visual  range  (extinction)  measurements  are 
summarized  in  Fig.  13.  It  shows  clearly  that  there  is  a 
significant  improvement  of  visibility  in  a~dense  smoke  (particle 
concentrations  surpassed  often  500,000  cm”3).  The  time  change  in 
light  ixtinction  is  measured  in  relative  units  (mV/min).  There 
is  not  much  difference  in  light  extinction  caused  by  coarse  fiber 
and  fine  (close  packed)  scavengers  what  supports  our  finding 
during  the  wind  tunnel  testing.  In  spite  of  a  large  spreading  of 
data  points  there  is  a  markable  difference  between  the  extinction 
background  (measured  before  the  scavengers  were  dropped),  marked 
by  a  dashed  line  in  Fig.  13,  and  the  measured  extinction  after 
the  scavengers  passed  the  sensitive  volume  of  the  beam  emitted 
from  a  light  source.  LED  columnated  source  (  X  •  6700  A)  formed  a 
light  beam  of  1.5  cm  diameter  and  80  cm  length.  The  light 
intensity  was  monitored  by  Bell  and  Howell,  Model  529,  Photodiode 
Detector  with  a  Transmission  Meter.  During  a  series  of  46 
measurements  the  most  significant  effect  (AV  >  10  mV/min)  was 
observed  immediately  after  the  first  scavenger  qroup  was  dropped 
and  lasted  approximately  for  1  hour.  The  number  of  scavengers 
dropped  usually  during  5  minute  Intervals  varied  between  10  and 
150  scavengers.  There  is  a  strong  indication  that  the  highest 
scavenger  dropping  rate  (30  scavenqers/min )  did  not  cause  the 
most  significant  change  in  light  extinction.  In  general, 
dropping  of  scavengers  in  smaller  groups  with  long  intermissions 
was  more  effective.  This  indicates  that  the  scavenger 
hydrodynamic  interactions  and  air  entrainment— which  lowers 
considerably  the  collection  efficiency  of  scavengers — might 
affect  considerably  the  experiments  in  a  smoke  chamber  of  limited 
volume.  For  this  reason  it  is  not  possible — at  this  time — to 
describe  by  a  theoretical  model  the  nature  of  smoke  particle 
deposition  on  coarse  or  fine  fiber  scavengers.  A  rough 
calculation,  based  on  the  formula  for  visual  range  in  a  foo 
(Podzimek,  1981) 
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with  e  *  0.02,  Kext  «  2.0,  n  (total  particle  concentration)  ■ 
500,000  cm-3  ancj  particle  size  distribution  np  (dp),  leads  to  L  of 
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the  order  of  10  cm.  This  is  consistent  with  the  observation  in 
the  smoke  chamber.  An  improvement  of  visual  range  by  a  factor  of 
1.3  or  1.4  seems  to  be  feasible  if  a  sufficient  number  of  hiqhly 
effective  scavengers  will  be  used.  For  instance,  from  the 
laboratory  studies  with  coarse  fiber  scavengers  (0.5”  x  0.5”)  the 
scavenger  surface  density  of  0.5  sc./cm^  would  mean  a  35.8% 
improvement  in  visual  range.  The  evaluation  of  electron 
micrographs  of  particles  deposited  on  coarse  and  fine  fiber 
scavengers  falling  in  the  smoke  chamber  reveals  that  the  fine 
fibers  collect  mainly  particles  smaller  than  0.5  %im,  whereas 
coarse  fibers  capture  also  micron  size  particles.  This  supports 
the  investigation  made  in  a  wind  tunnel  mentioned  earlier.  The 
dual  character  of  smoke  particle  deposition  on  coarse  fibers  is 
demonstrated  in  Fig.  12.  It  was  also  observed  that  TiCl.  smoke 
particles  deposited  on  the  obverse  side  of  the  fiber  around  the 
stagnation  point  (line)  form  chain-like  aggregates  probably  due 
to  the  combined  effect  of  electrostatic  and  inertial  force. 

The  experiments  in  a  cylindrical  smoke  chamber  showed  also 
the  importance  of  investigating  the  effect  of  hydrodynamic 
interaction  and  dispersion  of  a  large  guantity  of  scavengers 
released  simultaneously.  The  scavengers  dispersed  from  a  point 
form  a  cone  the  shape  of  which--besides  other  parameters--depends 
on  the  induced  air  circulation  in  a  specific  experimental  space. 
For  this  reason  comparative  measurements  with  the  same  type  of 
scavengers  have  been  made  simultaneously  in  the  smoke  chamber 
(0.0  mJ  volume,  1.6  m  high)  and  in  the  UMR  gymnasium.  There  a 
known  quantity  of  scavengers  was  dropped  at  a  height  of 
approximately  10  m  and  the  scavenger  position  on  the  floor  was 
evaluated  in  reference  to  the  center  line  of  the  "dispersion" 
cone. 


The  main  results  of  this  investigation  are  summarized  in 
Table  I,  where  the  following  parameters  are  calculated j  The  mean 
scavenger  velocity,  standard  deviation  of  deposited  scavengers 
from  the  cone  axis,  mean  "radius"  of  circle  within  which  were 
found  50%  of  the  total  number  of  scavengers  and  the  coefficient 
of  scavenger  dispersion.  One  can  observe  a  permanent  difference 
between  the  parameters  calculated  in  the  chamber  and  in  the 
gymnasium.  In  the  first  case  the  scavenger  settling  path  was 
approximately  1.6  m  and  in  the  second  one  10  m.  Settling 
velocities  for  the  same  type  of  scavengers  are  higher  in  the  tank 
than  in  the  gymnasium.  To  this  is  related  the  scavenger 
dispersion  on  the  floor  and  the  mean  "radius”  of  the  circle 
circumscribed  to  the  area  where  501  of  scavengers  are  deposited, 
?.  For  the  height  of  the  dropping  point  H,  for  the  coefficient 
of  scavenger  disnersion,  Dg,  and  scavenger  settling  velocity,  Vg, 
the  7  was  calculated  from  the  relationship, 

r  ■  (2Dst)1/2  (H  •  Ds)1/2  Vs’1/2  .  (11) 

Practically  all  data  indicate  larger  mean  radii,  r,  and 
dispersion  coefficients,  Dq ,  in  the  gymnasium.  The  explanation 
might  be  related  to  a  minimum  path  length  along  which  the  steady 


scavenger  motion  will  develop  and  the  strong  effect  of  induced 
air  velocities  in  the  smoke  chamber.  The  dispersion  coefficient 
was  found  to  be  very  large  in  cases  where  the  shape  of  scavengers 
supported  their  sliding  or  rotation  along  their  settling  path 
(e.g.j  paper  disks  and  rectangles  in  Figs.  14  and  15)  in 
comparison  to  steady  settling  of  mesh  type  scavengers  in  calm  air 
(Fig.  16).  Each  "Distance  Zone”  in  Figs.  14,  15  and  16  is  a 
multiple  of  27.4  cm. 

TABLE  I 
COMPARISON 

OF  SCAVENGER  DISPERSION  PARAMETERS 


SCAVEN. 

TYPE 

FALL  VELOCITY 
(cm/s) 

TANK  CYMN 

STAND. 

TANK 

DEV. 

GYMN 

50%  RADIUS 
(cm) 

TANK  GYMN 

COEFF.  DTFFnSTOf 

TANK  GYMN 

Coarse  rect. 
"dust  magn. " 

191.81 

163.66 

15.59 

4.79 

12.122 

43.293 

68.92 

159.6: 

Small  coarse 
rect.  "dust 
magnet" 

192.81 

166.35 

15.89 

4.79 

12.293 

33.790 

70.87 

98.8' 

Coarse  fine 
fiber  "dust 
magnet" 

176.35 

173.40 

8.79 

10.78 

33.000 

97.97 

Small  fine 
fiber  "dust 
magnet" 

146.05 

146.05 

146.05 

120.90 

18.29 

18.29 

18.29 

10.78 

19-308 

14.529 

15.52 

45.755 

133.13 

75.38 

129.1! 

Paper  disc 
r»0.6  cm 

150.76 

143.35 

9.40 

10.17 

33.840 

197.28 

516.10 

2892. 8( 

Paper  rectang. 
0.13x0.3  cm 

177.31 

177.31 

177.31 

175.70 

15.96 

15.96 

15.96 

11.95 

34.837 

36.542 

33.626 

130.565 

526.11 

578.87 

490.191 

1559. 3( 

Paper  square 
Paper  rectanq. 
0.13x0.3  cm 

177.31 

177.31 

177.31 

173.93 

175.70 

15.96 

15.96 

15.96 

14.976 

11.95 

34.837 

36.542 

33.626 

47.843 

130.565 

526.11 

578.87 

497.191 

207.3' 
1559. 3i 

Paner  square 

173.93 

14.976 

47.843 

207 . 3, 

0.11x0,14  cm 


Fig.  14  Number  of  deposited  circular  punch  card* 

(d  ■  0.66  cm)  in  dependence  on  the  distance 
from  the  vertical.  The  acavengera  were 
dropped  from  the  plateau  9.60  m  above  the 
floor. 
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Fig.  15  Number  of  deposited  paper  rectangles  (0.13  x 

0.3  cm)  in  dependence  on  the  distance  from  the 
vertical.  Each  circular  zone  width  on  the  floor 
is  27.4  cm.  The  scavengers  were  dropped  from 
the  plateau  9.60  m  above  the  floor. 


Fig.  16  Number  of  deposited  fine  fiber  plates  (1.25  cm 

x  1.25  cm  x  0.6  cm)  in  dependence  on  the  distance 
from  the  vertical.  Each  circular  zone  width  on 
the  floor  is  27.4  cm.  The  dropping  plateau  was 


9.60  m  above  the  floor. 


4. 


INVESTIGATION  OF  THE  FALLING  SCAVENGER  ZONE 


A  large  quantity  of  scavengers  induces  downdrafts  inside  a 
scavenger  zone.  These  downdrafts  and  related  horizontal  air 
motion  affect  the  scavenger  collection  efficiency  and  motion,  and 
potentially,  the  transport  of  smoke  particles  towards  the  ground. 

During  the  past  three  years  the  two-dimensional  model,  which 
in  the  first  stage  does  not  include  the  thermal  instability 
terms,  was  developed  (Podzimek,  1983?  Podzimek  and  Smith,  1986). 
The  model  is  suitable  for  investigation  of  the  impact  of 
different  scavenger  mixing  ratios  (concentrations),  isotropic 
turbulent  exchange  coefficient  and  release  height  on  the 
evolution  of  the  scavenger  cloud.  Recently  the  model  was 
extended  to  include  the  effect  of  slight  horizontal  wind  on 
scavenger  dispersion  and  induced  air  motion.  The  basic  set  of 
equations  with  the  boundary  and  initial  conditions  is  described 
in  Appendix  E. 

Summarizing  the  results  of  the  investigation  of  the 
evolution  of  the  falling  scavenger  zone  with  its  potential  impact 
on  the  smoke  particulate  transport  to  the  ground,  one  finds  the 
two-dimensional  model  suitable  for  describing  the  role  of  the 
main  parameters  affecting  the  zone  dynamics.  Apparently  the  most 
important  parameter  is  the  scavenger  mixing  ratio  (g  sc./g  air). 
For  non-spherical  scavengers  it  would  be  necessary  to  introduce 
an  "effective  mixing  ratio"  which  will  account  for  the 
aerodynamic  effect  of  non-spherical  collectors.  The  change  of 
mixing  ratio  from  0.01  g  sc./g  air  to  0.04  g  sc./g  air  is 
directly  proportional  to  the  speeding  up  of  the  scavenger 
dispersion  through  the  calm  atmosphere  (Fig.  17).  The  lowest 
position  of  the  zone,  z,  and  the  maximum  width,  x,  for  both 
mixing  ratios  are  plotted  in  Fig.  18  as  a  function  of  time.  The 
dispersing  cloud  generates  further  the  induced  motion  downwards 
(velocity  component,  w  in  Fig.  19)  and  entrains  the  air  from 
outside  with  the  horizontal  component,  u  (in  Fig.  20).  In  this 
case  scavengers  were  dropped  at  2  km  altitude  and  were  settling 
with  the  fall  velocity  of  1.5  m  s”  .  The  mo3t  intense  downdraft 
developed  4-fter  300  seconds  and  when  the  scavenger  zone  (of 
original  dimensions  400  m  x  400  m)  reached  the  ground  the  induced 
velocity  surpassed  13%  of  that  of  the  scavengers.  Very 
interesting  is  also  the  pattern  of  horizontal  velocity  isolines 
(Fig.  20)  which  demonstrates  the  intense  entrainment  of  the  air 
when  scavenger  zone  is  close  to  the  ground.  This  air  mass 
exchange  contributes  also  to  the  increasing  visibility  at  the 
ground.  The  transport  of  smoke  particles  downwards  is 
accompanied,  however,  by  the  decreasing  efficiency  of  inertial 
deposition  of  larger  smoke  particulates  on  the  scavengers, 
because  the  collection  efficiency  is  proportional  to  the  velocity 
difference  between  the  scavenger  and  aerosol. 

The  value  of  turbulent  exchange  coefficient  does  not  affect 
much  the  air  motion  around  and  inside  of  the  falling  scavenger 
zone.  Decisive  for  the  induced  air  motion  is  the  time  available 
for  the  full  evolution  of  the  air  circulation  and  the  initial 
shape  of  the  scavenger  zone.  The  following  combinations  of  the 
horizontal  and  vertical  dimensions  of  the  scavenger  zone  have 
been  used:  800  m  x  200  m?  500  m  x  320  m?  400  m  x  400  m  and  320  m 
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Fig.  18  The  mean  heights/  z,  of  the  falling  scavengers 
front  and  the  maximal  widths,  bmax  *  x,  as  u 
function  of  time  for  g  •  0.001  and  0.04  g  sc./g 
air.  Other  parameters  are  the  same  like  in  Fig. 
17. 
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19a, b,c  Impact  of  the  400  x  400  m  falling  scavenger 
zone  on  the  induced  w  -  velocity  components. 


U  /laid  (a/a)  V  •  1.3  a/a  •  200  aaa 

Fig.  20a 
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x  500  m.  The  most  significant  effect  have  thick  zones  (320  m  x 
500  m,  400  m  x  500  m)  which — for  scavenger  settling  velocity  of 
Vg  »  1.73  m  s~  --induced  downdrafts  surpassing  0.2  m  s”1  (Fig. 

21) .  Induced  downdrafts  by  thin  zones  (B00  m  x  200  m)  of 
scavengers  released  initially  at  an  altitude  of  1.5  km  induced 
velocities  at  the  ground  only  slightly  overriding  0.1  m  s*1  (Fig. 

22)  . 

The  numerical  modeling  of  the  evolution  of  the  falling 
scavenger  zone  at  slight  side  winds  concentrated  mainly  on 
"spherical"  scavengers  dispersed  at  a  high  altitude  (7  km). 

Their  settling  has  been  accompanied  by  induced  downdrafts  and 
considerable  deformation  of  the  zone  due  to  the  effect  of 
horizontal  velocity  (u  ■  1  m  s-i)  constant  in  the  whole  field. 

The  model  with  boundary  conditions  outlined  in  the  Appendix  E  can 
be  applied  to  other  velocity  profiles  and  to  the  changing  air 
density  with  altitude.  The  dramatic  evolution  of  the  scavenger 
cloud  after  a  long  time  of  settling  (Fig.  23)  indicates  the 
necessity  of  investigating  the  scavenger  dispersion  at  several 
typical  situations  in  the  troposphere  and  in  the  planetary 
boundary  layer  (e.g.#  stable  and  unstable  thermal  stratification) 
if  larger  areas  are  intended  to  be  covered  by  scavengers. 
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Fig.  22a, b,c  Impact  of  the  initial  shape  of  the  zone  (800  m  x 
200  m)  on  the  zone  evolution.  Calculated  for 
q  =  0.01  g  sc./g  air,  Hp  *  1  km,  V_  =  1.73  m/s, 
Kfc  =  1000  m2s.  K  & 

40 


m-*m t  m 


Fig.  23a, b,c  Falling  scavenger  zone  evolution  at  a  horizontal 
wind  1.0  m/s. 
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I  5.  DISCUSSION  OF  THE  MAIN  RESULTS 


The  main  goal  of  this  project  was  to  determine  the 
feasibility  and  ef fectiveneee  of  a  ecavenging  operation  for 
improving  visibility  in  a  military  smoke  cloud.  The  previous 
investigations  indicated  higher  scavenging  efficiency  of 
nonspherical  collectors  for  smoke  aerosols  with  particle  sizes 
from  0.2  ym  to  several  ym.  Different  shapes  of  scavengers 
(disks,  squares,  rectangles)  and  their  motion  do  not  have  a 
decisive  impact  on  collection  efficiency  which  rarely  surpassed 
1%  for  a  disk  and  Re  up  to  20  according  to  numerical  models 
(Podzimek,  1981).  Experimental  data  were  almost  one  order  of 
magnitude  smaller  than  the  theoretical  values  for  smoke  particle 
radii  around  0.3  ym.  Collection  efficiencies  of  electrically 
charged  disks  varied  between  3%  and  101,  however,  the  problem  of 
how  to  maintain  the  high  electric  charge  on  collectors  during 
small  scale  experiments  remained  unresolved.  Finally,  it  was 
found  that  for  smoke  particle  sizes  within  the  "Greenfield  gap" 
the  most  effective  scavengers  are  mesh  or  grid  type  models, 
mainly  if  "self-charging"  electrete  type  materials  are  used. 
Concerning  the  hydrodynamic  interaction  of  a  population  of 
scavengers  falling  simultaneously  in  the  "Falling  Scavenger  Zone 
was  recommended  to  follow  closely  the  zone  evolution  in  the 
future  with  regard  to  the  results  included.  In  this  report  wo 
would  like  to  call  attention  to  the  new  findings  made  during  the 
past  three  years  and  to  the  results  obtained  while  using  improved 
equipment  or  numerical  models.  Also,  we  would  like  to  address 
the  queetion  of  practical  applicability  of  our  study  for  medium 
or  large  scale  experiments  in  the  field. 

Collection  efficiency  of  thin  oblate  spheroids  was  extended 
to  Re  •  80.  At  Re  higher  than  80  our  numerical  model — which  is 
based  on  creeping  flow  equations— failed  to  yield  reasonable 
results.  Still  remains  to  explain  the  shape  of  the  collection 
efficiency  curve  for  0.1  H  Re  H  5.0.  In  this  transitional  region 
starts  to  develop  the  instability  on  the  reverse  side  of  a  disk 
which  has  an  impact  on  particle  deposition  through  phoretical 
forces.  The  additivity  of  the. di f fusion  and  phoretic  forces  was 
questioned  by  Carstens  and  Martin  (1982),  however,  it  can  be 
assumed  that  the  effect  of  these  forces  will  be  smaller  in  a 
smoke  cloud  than  in  water  droplet  or  ir.  a  mixed  cloud  (Podzimek, 
1984) . 

The  study  of  particulate  deposition  on  cylinders  was 
intended  to  confront  some  of  our  results  of  wind  tunnel 
measurements  with  a  simple  theory  of  deposition  on  a  single  fiber 
or  a  pair  of  hydrodynamically  interacting  fibers.  Compared  to 
the  known  solutions  of  particulate  deposition  in  the  potential 
flow  (Albrecht  1931;  Sell  1931;  Glauert,  1940;  Langmuir  and 
Blodget,  1946;  Davies  and  Peetz,  1955)  the  presented  solution 
might  be  more  accurate  and  very  close  to  the  results  of  the 
Davies  and  Peetz  study.  It  can  also  be  used  for  the  study  of 
particle  deposition  on  two  cylinders  close  together  and  yield 
quite  satisfactory  results  for  large  Stokes  numbers  of 
particulates  (sizes  larger  than  1  ym  and  densities  ^  2.0  g 
cm-3).  For  these  conditions  the  particle  deposition  rate  might 
be  close  to  the  deposition  in  viscous  flow  solved  in  a  numerical 
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model  by  Schlamp  (1977).  Hie  model  deduced  for  the  deposition  of 
cloud  droplets  on  columnar  ice  crystals  was  available  to  our 
group.  It  was  applied  to  the  deposition  of  smoke  particulates 
smaller  than  cloud  droplets  in  our  group  with  mixed  success. 

After  we  became  familiar  with  the  Brown's  multi-fiber  model 
(Brown,  1983)  which  is  described  in  Appendix  D  we  tried  to  use  it 
for  particle  deposition  on  parallel  or  staggered  fiber  array.  In 
our  opinion,  the  model  has  to  be  improved  in  order  to  describe  in 
a  more  realistic  way  the  conditions  inside  a  filter  type 
scavenger.  On  the  other  side,  it  has  great  potential  if  applied 
to  our  grid  type  scavengers  or  thin  not  very  dense  packed  filter 
type  scavengers. 

This  study  supports  the  conclusion  from  the  previous  report 
that  the  effect  of  an  electric  charge  is  significant.  The  only 
Coulomb  force  interaction  is,  however,  an  oversimplification  of 
the  whole  process.  The  new  calculation  supports  also  the 
previous  conclusion  that  the  collection  efficiency  is  lower  for 
particle  radii  larger  than  5.0  urn  if  compared  to  the  uncharged 
collector.  We  might  speculate  about  the  neglected  effect  of 
image  forces.  Finally,  the  discrepancy  between  the  calculated 
and  measured  values  of  collision  efficiency  (Fig.  2)  might  be 
explained  by  the  fact  that  the  model  assumes  a  single  first 
collision  of  the  smoke  particle.  The  second  and  all  following 
particles  are  impacting  the  collector  as  if  the  preceeding 
collisions  would  not  happen  or  if  the  aerodynamic  and 
electrostatic  field  around  the  collector  would  remain  the  same. 
Conclusions  similar  to  those  mode  for  thin  oblate  spheroid 
bearing  an  electrostatic  charge  can  be  made  for  a  cylinder 
(Podzimek  and  Martin,  1984). 

Comparing  the  theoretical  results  of  particle  deposition  of 
collectors  of  simple  geometry  (disks,  cylinders)  with  the  wind 
tunnel  measurements  or  with  the  observed  deposition  of 
particulates  on  scavengers  falling  in  a  smoke  chamber  we  observe 
usually  smaller  deposition  on  exposed  scavengers.  The 
explanation  might  be  due  to  the  imperfect  modeling  of  particulate 
behavior  in  the  boundary  layer  (particle  rotation  in  a  velocity 
gradient,  wake  effect).  Considerable  part  of  particulates  might 
rebounce  from  the  collectors  surface.  The  last  hypothesis  is 
supported  by  the  fact  that  charged  collectors  showed  a  deposition 
comparable  with  theoretical  results.  The  main  discrepancy  might 
be,  however,  explained  by  the  used  technique  of  deposited 
particle  counting  on  the  collectors  surface.  Usually  the  smoke 
particles  are  piled  up  or  make  some  larger  aggregates--which  in 
principle  will  increase  the  collection  ef f iciency--but  it  is  very 
difficult  to  count  all  particles  in  the  electron  micrographs. 

Many  of  the  small  particles  will  not  be  seen  in  the  proximity  or 
under  the  large  ones.  Several  of  the  trials  to  check  the 
reproducibility  of  the  particle  counting  led  to  the  conclusion 
that  the  mean  errors  of  counting  might  override  +  30%  of  the 
counted  value. 

One  of  the  most  significant  results  of  the  experiments  seems 
to  be  the  high  smoke  particle  deposition  on  the  reverse  side  of 
scavengers.  In  the  case  of  a  disk  it  might  be  comparable  to  the 
deposition  on  the  obverse  side  which,  unfortunately,  cannot  be 
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described  by  numerical  models  at  Re  >  10  (for  a  disk).  The 
scavenger  motion  or  oscillation  will  not  effect  significantly  the 
deposition  of  particulates  on  scavengers  considered  in  this 
report . 

The  experiments  in  cylindrical  smoke  chamber  were  marked  by 
a  large  scatter  of  data  from  the  measurement  of  the  light 
extinction  after  the  scavengers  were  dropped  in  a  dense  TiCl. 
smoke.  The  arrangement  of  the  experiment  was  the  same  like  that 
described  in  our  previous  report  (Podzimek,  1981).  More  emphasis 
was  put.  however,  on  the  reproducibility  of  the  measurements. 

The  generation  of  TiCl.  aerosol  was  controlled  within 
approximately  +  20%  of  deviation  of  the  individual  points  on  the 
size  distribution  curve  during  the  total  duration  of  the 
experiment  (usually  more  than  1  hour).  The  large  fluctuation  of 
signals  in  the  past  was  explained  by  the  different  mode  of 
scavenger  dropping  (frequency  and  initial  position  of  scavengers) 
and  the  low  probability  of  hitting  the  sensitive  volume  of  the 
light  beam.  Further  was  observed  that  large  quantity  of 
scavengers  dropped  at  the  same  time  did  not  cause  a  markable 
change  in  light  extinction.  This  is  explained  by  the  induced 
downdraft  in  the  relatively  small  smoke  chamber.  It  is 
hypothesized  that  the  downdraft  is  reducing  the  velocity 
difference  between  smoke  particle  and  scavenger,  and  in  this  way, 
also  the  collection  efficiency  of  scavengers  will  be  smaller. 

The  observation  of  scavenger  interaction  and  formation  of  a 
dispersing  cone  in  the  smoke  chamber  1.6  m  high  led  to  the 
question  of  how  to  transfer  the  data  from  the  chamber  into  the 
real  atmosphere.  Experiments  in  the  UMR  gymnasium  showed  clearly 
that  .the  difference  in  the  dispersion  angle  depends  on  scavenger 
type,  size  (mass)  and  motion,  but  also  on  tho  rate  of  scavenger 
dropping.  Scavengers  required  a  certain  time  (several  seconds) 
for  forming  a  steady  dispersing  cone  in  spite  of  the  fact  that 
their  relaxation  time  for  resuming  a  steady  motion  (of  a  single 
scavenger)  is  of  the  order  of  tenths  of  a  second.  We  believe, 
that  the  observed  dispersion  cones  or  coefficients  of  scavenger 
dispersion  in  the  calm  atmosphere  are  useful  parameters  for 
practical  use  in  the  field.  The  data  were  also  used  for 
estimating  the  overall  effect  of  SPADES  technique  for  clearing  a 
military  smoke  cloud.  A  logical  continuation  of  this  study  ought 
to  be  the  experiments  with  the  most  effective  (electrete  grid  or 
filter  type)  scavengers  at  a  slight  side  wind  in  the  atmosphere 
or  in  an  aerodynamic  wind  tunnel  of  large  cross  section. 

The  methodology  of  a  successful  application  of  scavenger 
technique  for  clearing  of  a  smoke  cloud  has  to  be  combined  with 
an  assessment  of  the  clearing  effect  through  induced  downdrafts, 
with  regard  to  the  dispersion  of  scavengers  from  a  helicopter, 
airplane,  or  parachute  the  2-D  model  simulating  the  effect  of  a 
"line  source"  seems  to  be  an  appropriate  approximation.  Its 
extension  for  simulating  the  scavenger  dispersion  at  a  slight 
side  wind  and  for  atmosphere  with  variable  air  density  has. been 
already  done.  Especially  the  study  of  the  effect  of  1  ms'1 
horizontal  wind  is  an  important  contribution  to  the  study  of 
clearing  polluted  atmospheric  boundary  layer  by  artificial  or 
natural  (rain  drops,  ice  crystals)  scavengers.  The  further 
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evolution  of  this  model  could  go  via  inclusion  of  the  thermal 
stratification  of  the  atmosphere  {similar  study  with  different 
goals  has  been  made  by  Girard  and  List,  1975)  and  application  of 
a  specific  wind  velocity  profile  in  the  boundary  layer. 

Scavenger  mixing  ratio  seems  to  be  one  of  the  most  important 
parameters  directly  related  to  the  generated  downdrafts.  This 
leads  to  a  logical  conclusion  that  for  an  effective  clearing  of 
smoke  cloud  have  to  be  combined}  That  collecting  the  particles 
on  special  type  of  scavengers  and  that  supporting  their  transport 
towards  the  ground.  The  second  technique  would  not  require  the 
use  of  any  type  of  effective  scavengers  and  the  positive  effect 
could  be  obtained  by  dispersing  dry  ice  grains,  water  droplets, 
or  provoking  artificial  rain  over  the  area  at  favorable 
meteorological  situation.  Other  important  parameters  for 
inducing  an  effective  downdraft  are  the  settling  velocity  of 
scavengers,  the  minimum  time  required  for  downdraft  evolution  and 
the  initial  thickness  of  the  scavenger  zone.  They  can  be  easily 
controlled  during  medium  or  large  scale  experiments. 

In  compliance  with  point  4  of  part  "8.  Identification  of 
the  Main  Goals"  in  our  proposal  to  the  program  DAAK-11-83 — -which 
we  intended  "to  establish  a  solid  basis  for  the  application  of 
scavenging  technique  in  the  field..."  we  tried  to  quantify  the 
impact  of  the  main  factors  contributing  to  the  evolution  of 
downdrafts  in  the  falling  scavenger  zone.  These  factors  depend 
also  on  the  applied  technology  of  scavenger  dispersion:  Several 
scenarios  were  suggested  for  electrete  fiber  dispersion  from  an 
airplane,  helicopter  or  a  parachute.  Three  groups  of  four  senior 
students  studied  (without  knowing  the  details  about  the  program 
DAAK-fll-83)  the  technology  of  scavenger  release  and  dispersion. 
The  reports— available  in  our  office  include  the  details  on  the 
scavenger  generators,  estimates  of  generators  output  and  the 
amount  of  scavengers  necessary  for  30%  improvement  of  visibility 
at  the  ground.  These  numbers  (in  Table  II)  are  based  on  the 
measured  data  in  the  wind  tunnel  and  in  the  smoke  chamber  and  do 
not  include  the  effect  of  induced  downdrafts. 

Particle  removal  efficiencies  in  Table  II  represent  the  most 
conservative  numbers  which  were  deduced  from  a  model  (Podzimek, 
1979)  and  wind  tunnel  measurements  with  single  fibers.  The 
effect  of  fiber  crossing  in  the  mesh— which  might  enhance  the 
collection  ef f iciency— is  not  considered. 

The  total  effect  would  be  certainly  more  favorable  depending 
on  the  choice  of  scavengers  for  the  most  effective  downdraft. 

Heavy  scavengers  with  compact  structure  will  be  more  suitable  for 
a  thin  layer  of  smoke  in  reference  to  downdraft  generation. 
Apparently  the  combination  of  different  techniques  during  one 
action  will  be  the  most  appropriate  way  for  clearing  a  military 
3moke  cloud  in  our  opinion. 
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TABLE  II 


TOTAL  WEIGHT  OF  SCAVENGERS  NECESSARY  FOR  IMPROVING 
THE  VISUAL  RANGE  OF  30%  OVER  AN  AREA  OF  100  m  X 
100  m  IN  Ti02  SMOKE 


Scavenger  Type 
size  [cm] 


Particle  Removal  Efficiency 
Scav.  Weight  in  %  for  particle  radii  [jwn] 
[gl  <0.25  0.50  0.75  1.00  1.25 


Weight  of 
Scavengers 
(kg) 


I  Smail 

coarse  fiber 

|  1 . 3  x 

1.3  x  0.71.4 

0.023 

4.7 

3.2 

1.1 

0.9 

0.3 

815.0 

i  Large 

3.1  x 

i 

coarse  fiber 

3.1  x  0.024 

0.125 

4.7 

3.2 

1.1 

0.9 

0.3 

1,632.0 

.5  3.2  0.6 


5  3.2 


600.0 
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6.  CONCLUSION 


The  program  of  three  year  effort  was  formulated  in  our 
proposal  in  paragraph  “7.  Time  Schedule"  and  in  "8. 
Identification  of  the  Main  Goals".  All  points  of  our  program 
were  followed  in  accordance  with  our  proposal,  however,  with 
different  intensity.  For  example,  interaction  of  diffusional  and 
phoretical  forces  was  not  followed  in  much  detail  when  it  was 
clear  that  diffusion  and  dif fusiophoresis  will  not  play  a  major 
role  in  smoke  particle  (0.1  <  r  <  5.0  ym)  deposition  on 
insoluble  fibers.  Also  the  effect  of  oscillating  scavengers  had 
an  orientational  character  due  to  the  little  effect  of  the 
observed  frequencies  and  scavenger  motion  amplitudes  on  the 
collection  efficiency  of  planar  scavengers. 

Most  effective  were  the  grid  type  scavengers  ana  filter  type 
scavengers,  especially  those  made  of  electrete  fibers.  The 
physical  mechanism  of  their  action  is  not  completely  known, 
however,  their  collection  efficiencies  (referred  to  the  ideal 
cylinder  type  fiber)  will  surpass  10%.  This  finding,  based 
mainly  on  the  wind  tunnel  experiments,  seems  to  be  supported  also 
by  the  results  from  the  smoke  chamber.  In  the  future  attention 
should  be  paid  to  the  appropriate  combination  of  fibers  with 
different  geometry  and  composition  in  order  to  reach  a  maximal 
effect  in  polydisperse  smoke  particle  removal. 

For  medium  and  large  scale  experiments  in  the  field  the 
knowledge  of  dispersion  characteristics  of  a  large  population  of 
falling  scavengers  seems  to  be  of  utmost  importance.  Experiments 
in  large  scale  (10  m  settling  path)  showed  clearly  that  each  type 
of  scavengers  needs  a  specific  time  for  developing  the  "typical" 
shape  of  a  dispersion  cone  and  that  the  experiments  in  small 
smoke  chambers  must  be  corrected. 

During  the  study  of  direct  scavenging  of  particulates  was 
discovered  that  the  air  downdraft  induced  in  the  falling 
scavenger  zone  might  be  another  factor  contributing  to  the 
effective  transport  of  smoke  particles  to  the  ground.  The 
magnitude  of  the  downdraft  depends  on  the  scavenger  mixing  ratio 
and  settling  velocity.  The  time  required  for  the  fully  developed 
downdraft  (surpassing  often  10%  of  scavenger  settling  velocity) 
and  the  initial  thickness  of  scavenger  zone  are  other  important 
parameters  determining  the  evolution  of  induced  air  circulation 
and  transport  of  smoke  particles  to  the  ground.  The  appropriate 
combination  of  direct  scavenging  with  the  induced  downdraft 
technique--which  can  be  accomplished  by  simple  water  droplet 
dispersion  or  cloud  seeding--might  be  an  effective  technique 
especially  for  smoke  clouds  formed  during  very  stable  thermal 
stratification  in  the  atmospheric  boundary  layer. 
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APPENDIX  A 


CALCULATION  OF  PARTICLE  TRAJECTORIES  FOR  POTENTIAL  FLOW 


For  sufficiently  high  Rc  the  model  by  Bergrun  (1951) 
(potential  flow  around  a  cylinder)  and  for  small  Re  the  model  by 
Schlamp  (viscous  flow)  was  applied.  The  potential  flow  assumes 
that  the  inertial  forces  outweigh  the  frictional  forces  in 
particle  deposition  and— if  the  effect  of  boundary  layer  can  bo 
neglected — that  a  good  description  of  the  particle  deposition 
will  be  obtained  for  the  front  side  of  the  cylinder.  The  actual 
flow  on  the  reverse  side,  however,  deviates  much  from  reality. 


The  basic  equations 
of  the  ideal  flow  pace  a 


for  the  radial  and  tangential  components 
cylinder  of  radius  Rc  are  (Liu,  1935) 


vr  “  u»  ll-4)  cos  6  5  ?  ‘If1 
R  2 

ve  *  u.  Iin  1  =  ■  (ff> 


(i-i) 

(1-2) 


with  the  corresponding  stream  function 


tp=ursin0  (1 
00 


(1-3) 


From  the  definition  and  chain  rule  follows 


Ug  *  rdQ/dt  *  r (30/3 x ) ( 3 x/3 t ) +r (30/3y ) (3y/3 t )  ;  (1-4) 

ur  =  dr/dt  »  (3r/3x) ( 3 x/3 t ) + ( 3 r/3 y ) ( 3y/3 t )  .  (1-5) 

From  the  trigonometric  relation 
tan8  =  y/x  ;  9  =  tan-1(y/x) 

and 

r  =  (x2+y2)0,5  , 
can  be  deduced 
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3  6/3 


x  «  -yx“2/(l+(y/x)2)  *  -y/ ( x2+y 2 )  ; 


36/3y  *  x“l/(  l+(y/x) 2  ■  x/(x2+y2)  ; 

3  r/3  x  *  0.5(x2+y2)"0*52x»x(x2+y2)“0*5  ? 

3r/3y  ■  0.5(x2+y2)”°*^2y»y(x2+y2)"0,5  . 

Substituting  these  equations  into  equation  (1-4)  and  1-5),  one 
obtains 


uQ  »  -y(x2+y2)"0, 5vx+x(x2+y2)"0* 5vy 

-u  fl+R  2/(x2+y2) ly/(x2+y2)“°* 5  ; 


ur  *  x(x2+y2)-0,5vx+y(x2+y2)'°‘5vy 

*  uao(l-Rc2/(x2+y2)  ]  x/(x2+y2)-0*  5  ; 
-yux+xuy*-uao(l+Rc2/(x2+y2)]y  ? 

*ux+/VuJ1-Rc2/(*2+>'2,lx  * 

From  the  last  two  equations  were  calculated  velocities,  ux  and  uv 

(1-6) 

(1-7) 


u  »u  ( 1-R  2 ( x2-y2 )/( x2+y2 ) 2 ]  ; 

X  oo  C 


uy*-2UaoRc2xy/(x2+y2)2 


From  relation  y  =  rsin6  one  can  obtain  stream  function  easily  in 
dimensionless  way  from  equation  (1-3) 


H^yll-R  _ -/  (  x  2  +y  2  )  )  ; 


(  1-C  ) 


equations  (1-6),  (1-7),  (1-8)  can  be  nond imens iona 1 i zed  by 

setting 
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x'»x/Rc  ;  y'ay/Rc  ;  Ax'aAx/Rc  ;  Ay'aAy/Rc  f 
t'“tu«/Rc  :  u'-Ujj/u.  ;  u'-Uy/u^  ;  * '  uaoRc )  . 
with  these,  ju-»*:icns  (1-6),  (1-7),  (1-8)  become 


••1-.;  y ‘ 2 )/(x ' 2+y ' 2 ) 2  , 

(1-9) 

'  a-2x ' y ' / ( x ' 2+y ' 2  )  2  , 

( 1-10) 

■y ' ( 1-1/ ( x ' 2+y ' 2 ) 1  . 

(1-11  ) 

For  potential  flow  with  electric  force  the  particle 
trajectory  can  be  deduced  from  the  equation  of  particle  motion 
(velocity  components  vx,  v  )  in  which  the  gravitational  force 
will  be  neglected:  y 


dv/dt  *  -<CDRe/24K)(v-u)-Ceer  ,  (1-12) 

g*^r'2),  where  K*2pr  2u00/9vRc  is  the  Stokes  number, 
0,5/R  Ln  is  a  Coulomb  force  number,  i.e.  Cgor  is 
ess  Coulomb  force  between  two  point  charges  QA  and 

Qb  separated  by  a  distance  r  ■  r*Rc. 

The  transformation  of  Coulomb  force  number  and  electrostatic 
force  term  is  as  follows: 


C.-V  C/(u»2k 

aRd  f  S  oA0BK 
the  dlmensTortl 


Fe  = 


K0OAOQ/r‘ 


m. 


4rrp  P p/3 


hence 


V”?  ■  3K0QA00/(4»rp3pp3l 


( 1-13) 


g 

1  Coulomb  =  3  x  10 


1  OAl.fh  AO 


2  1 

j  »  1  dyne  *  cm  /(s.s.u.)*’  , 

let  Q  '  =  Q. / ( LP.  ) ,  Oo'  =  On/r_2,  the  electrostatic  force  term 

•  m  M  C  L5  D 

becomes 
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-Aryl  r./1 


I3W  )/(‘l'Trpppr2)-(Rc2L/u.2)e1.  . 

(30A’0BV(4ir)n2uiii/(0MOpRc)l0l5/[0prp2V,S«Rc|l0-5L/lu.2l"',0r- 

Let 

Vt3<VV/(4lT)l  l2u<jo/(9u0pRc)l0-5  * 

-OA'Q0'K°'5/(LRcmp)  (1-14) 

be  the  Coulomb  force  number.  Then  the  corresponding 
electrostatic  force  term  is 

Ceec  -  (2Cprp2u„/(9PRc)]-0'5-FEL/(u_2r'2)er  . 

At  this  point  Bergrun's  scheme  is  modified  to  get  the 
equation  of  trajectory  of  particles 

vxdvx/dx  «  -(CDRQ/24K) (vx-ux)-Ce  cos8  ,  (1-15) 

vydvy/dy  -  (CDRe/24K) (vy-uy)+Ce  sin6  ,  (1-16) 

dx/dt  *  _vx“(“e  cos®  <  dy/dt  *  -vy+C@  3in  , 

’V^Re’3  •  ''W2  *  (Vuy’2  •  <  1-17) 

The  velocity  and  position  of  a  particle  at  the  end  of  any 
interval  of  time  may  be  expressed  by  the  following  equations: 

vx ( n  +  1 ) =vx ( n ) ” ^ K ) ( vx ^ n j -ux ( n j ) At -Cecos 9  At  ,  (1-1C) 

Vy  (  n  +  1 )  *vy  (  n  )  "^  cDRc//^4K  ^  *  vy  (  n  )”uy  (  n  )  ^  ^t+Ces  in  9  At  ,  (1-13) 

VX(n+l )  =vx  (  n  )  t_3*5(CDRc//24K)  (vx(n)“ux(n)  5  t  2-0 . 5Cecos9At"  ,  (1-20) 

y(n  +  l)=vy(n)  t-0 . 5  ( CDRc/24 1< )  ( vy  (  n  )  -uy  (  n  )  )  1 2  +  0  .  SC^sinO  At2  .  (1-21) 
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Assume  the  initial  velocity  of  the  particle  to  be  the  same  as 
that  of  the  air  at  a  fairly  large  distance  forward  of  the 
cylinder;  then  one  can  estimate  the  second  velocity  of  the 
particle  both  in  the  X  and  Y  directions.  The  relevant  relation 
is  given  in  Appendix  B. 

vxl*°‘5{ (uxo"AxH)+l(uxO"AxH)2+4(uxlAyH+AxCecos9),°’5^  *  (1”22) 

vyl»0.5{ (uy0-AyH)+((uy0-AyH)2+4(uylAyH+AyCesine)]0,5}  .  (1-23) 

A  listing  of  the  computer  program  is  attached. 
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APPENDIX 

LISTING  OK  PROGRAM 


c  THIS  PROGRAM  COMPUTES  DROP  TRAJECTORIES  PAST  A  CYLINDER  IN 
c  POTENTIAL  FLOW  WITH  ELECTRIC  CHARGE  UN  CARTESIAN  COORDINATED 
c  RELAXATION  TIME:  TAU=(2’'RP*RP*RHOP)/ l 9-MEUA) 
c  STREM ( I )»U INF*ETA( I )*( I -RC'KC/ iPStU )  -*2+KTA( I )**2 ) ) 
c  STORE’S  NO  :  ST=TAU*UINF/RC 

c  LENGTH  UNIT  :  CM  -  BUT  IN  DIMENSIONLESS  WAY 

c  CHARGE  UNIT  :  ELECTROSTATIC  UNIT 


e 


REAL  MEl'.MASSP 

COMMON  MEU .  RHOA .  RHOP  .RC.RP.N.VXUOO)  ,  VY  (  -00  )  .  K  .  AK .  QAA ,  QBB ,  KENOC 
DIMENSION  VPSI(iOO) ,YETA(400) .PS  1(400  i  .F.TAUOO) ,RF.N0L(400) , 

5  C0EF(400)  .ELEF(400)  .RSQUOO)  ,I)IST(400)  .C0ST(400)  . 

S  SI NT l 400) 

CALL  INITT1960) 

CALL  TERN (3, 1024 D 
CALL  CHRSi::(4) 

CALL  DW I NDO (-6. 0,2. 0,-2. ,4. ) 

CALL  TWINDO( 100.900. 100.700) 

CALL  ANMODE 

READ l 5 . * ) RHOA , RHOP . nEU . N 
READ ( 5 . * ) RENOC . DC . DP . AR . QAA . QBB 
RP-0 . 5*DP 
RC-0 . 5 ••'■DC 


U I  SF*RENOC*MEU/  (  DC-'-RHOA ) 

WR ITE (6, 4) DC, DP, U INF, RHOP , AR , RENOC 
4  FORMAT ( IX .  ••'•■'DC» ’  .E12.5, IX, ’**DP*’  .E12. 
SRHOP*' ,F5 . 2 , IX , ,<nVRATIO  AR  *' ,£12.5. IX. 
WRITE (6 , 22) QAA, QBB 

22  FORMAT! IX, 'QAA3' . F6 . 2 , IX . ' QBB=’ ,F6.2) 
CALL  CIRCLE 


3  .  IX. ' ,>*UINF»'  . FI 0.4/ IX.  ’ 
’•'’•'RENOC  ■’ .F10.3) 


DO  55  K*l,  12 
CALL  ANMODE 
READ(  5  DFS 


WRITE  t  o ,80 )FS 

SO  FORMAT i  /  /  IX .  ’  ***FS«  ’  .  E 1 2 .  o  ) 
CALL  DISPlUISF. PSI.ETA.PS. ID 
CALL  MOYABSt 100. 300 D 


55  CONTINUE 


CALL  FINlTTlO . -000) 

STOP 

END 


ROUTINE  FOR  CALCULATING  T11F.  VELOCITY  OF  I'kLL  STREAM 
SUBROUTINE  FREST ( LINK . ! . PS ) 

REAL  MLU 

COMMON  MEU  ,  RHOA .  RHOP .  RC  .  KP  ,  \  X(  -‘JO  D  ,YY<-0<> )  ,  K  ,  AR  ,  QAA  ,  Ql) B  ,  RLNOC 
X--0 . I 


P 1  =  1  . 
p;=-fs 
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YY*0. 10286’'K 
CALL  MOVEA(-o.O.YY) 

DO  50  1*1  ,N 
X*X+0 . 1 
P3»X*X*1. 

P4=-FS*X*X 

CALL  SOLVE (PI , P2,P3,P4,Yl,Y2,Y3) 

CALL  DRAVA(X.Yl) 

A*X*X+Y  l'-Y  1 
D«X*X-Y1*Y1 
B»A*A 

VX( I )*1 . -D/B 
VY(I)*-2.»X*Yl/B 
50  CONTINUE 

CALL  M0VEAC-6.0.YY) 

RETURN 

END 

q  •Srin'rfrti  V# Vr  VnV  -ft  t  <r  VnV 
c  *>  •.WrfifrtWHhWrtKfrtV/ 

c.  .  . .ROUTINE  TOR  SOLVING  CORRESPONDING  Y  COORD  FOR  EACH  X  COORD 
SUBROUTINE  SOLVE! Pi .P2.P3 ,P4,Y1 ,Y2 ,Y3) 
F(YY)*P1*YY**3*P2*YY’'YY+P3*YY+P4 
G ( YY )*3 . *P I*YY*YY+2 . *P2*YY+P3 
EPS*1 .  E~0A 
1 I*-2 
JJ«0 
YI*0. 

15  Y=Yl+5. 

2  YI*Y-F (Y)yG(Y) 

JJ*JJ+1 

IF(JJ.EQ.  100)  GO  TO  10 
IF(ABS(F(YI ) ) -EPS)  4,3,3 

3  Y=YI 

GO  TO  2 

4  II=II*l 

iran  7,s.o 

1  Yl»Yl 
GO  TO  5 

5  Y2*YI 
GO  TO  5 

«  Y 3=YI 

5  IFl II .NE.  1  )  GO  TO  15 
GO  TO  12 

10  WRITE (6, 11  * 

11  FORMAT (.  IX,  ‘ THE  SOLUTION  OF  Y  CAN  NOT  CONVERGE'  ) 

12  CONTINUE 
RETURN 
END 

<;;;************* 

£  *.’?  V.**V  Vr*/r  *.V  Vf  V  V/V/'VV:V/Vr 

c. . . .ROUTINE  FOR  CALCULATING  THE  DISPLACEMENT  OF  PARTICLE 
SUB  ROUT I NE  D I SP ( U I  NT . PS  I , ETA . FS  ,  I ) 

REAL  MEU.NASSP 
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n  n  n  on  on  ooo 


COMMON  MEL* . RHOA . RHOP . RC . RP . N , VX ( 400 ) . VY (400 ) . K . AR .QAA.QBB, RENOC 
DIMENSION  VPSl(iOO)  .VETAUOO)  .PSl^OO)  .ETAUOO) . 

S  RENOL(400) ,C0EF(400) ,ELEF(400) ,RSQ(400) , 

$  DIST(.400) .COST! 400)  ,Sl.\Tl400) 

DELPS-0. 1 

psun^-b.o 

ETA(I)»0.  lOCSo'-K 
PSIO-PSH  n+DEI.rS 
CALL  FREST(.UINF,  I  ,FS) 

CALL  MOVEA  (  FS  I(.  1 ) .  ETA  Cl ) ) 

VPSlll)*VXll) 

VETA( 1 )*VY( 1 ) 

DEL£T=VY ( 1 1 /VX( 1 )*DELr S 
ETA(2)*£TA( 1 )+DELET 
Pl=3  Ul5‘i2o 

MASSP*4  .  /  3  .  *PI*RP**3*RHOP 
WRITE(6 , 90)MASSP 
90  FORMAT! IX, ’MASSP*’ ,E12.6> 

c _ CALCULATE  THE  REYNOLD  NUMBER  OF  FREE  STREAM 

R  f  NO P a R I ! 0  T I  \ F* R r*  2 .  /MEU 
VRITEco.  5  1DELET 
3  FORMAT U  X , ’ DE LET® ’  ,E12.o) 

WRITE! 6 , 8 1REN0P 
S  FORMAT !l X. 'RENOP-' ,FS. 3) 

c _ CALCULATE  THE  RELAXATION  TIME 

TAl  *2  .  »RP*RP«RI10P/  (9 . ’'MEU) 

ST»TAU*UINF/RC 
WRITER,  IDST.TAU 

11  FORMATUX,  'ST*'  ,EI2.6,  IX,  *TAU=’  ,E12  6) 

RES0L(1)*0. 

COEFU)*I. 

RSQU  )3PSI  ( I  2+ ETA (.  t  >**2 
DIST(1)»SQRT!RSQ(D) 

THE  X  AND  Y  COMPONENT  RELATED  TO  THE  ELECTRIC  CHARGE 


COST  1 1  )*rSIU)/DISTU  1 
SINT! 1 )*ETA( l )  /DISTl 11 

.FIND  THE  LENGTH  OF  CYLINDER 
CLL.NCsRC/AR 


. . .FIND  THE  MASS  OF  CYLINDER 
CMASS»3 .  14l5‘)2o*RC*RC*CLESG 
RMASS*MASSP/CMASS 
WRITE (6, 911 CLUNG, CMASS.RMASS 
FORMAT ( IX,  ’  •CLUNG* '  ,E  1 2  .  5  ,  IX  ,  '  CM  ’  .  IX 
5  '  •••"••  RMASS  = '  ,K12  5) 

IF  (KMAS5 *0.001  )  -»4,44,55 


.::;CMASS= '  .F.12. 5  ,  IX,  'GM'  / IX  , 


ELECTRIC  CHARGE  DATA  IN  DIMENSIONLESS  WAY 


44  OA-RC'C  LENG  QAA 
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QB*RP*RP*QBB 
VRITEio.  11  DQA.QB 

111  F0RMAT(1X, 'QA*' , £12.5, IX, ’QB-' ,E12.S) 

C 

c .  .  .  .  INITIAL  VALUE  OF  COULOMB  FORCE  NUNBER  AND  ELECTRIC  TERM 
c 

FE l*QA,'-'QB’;*STrt^0 . 5/  (RC*NASSP*CLENG) 

ELEF(  1  )=FF.I*C!.ENG/ (L'INF'**2*ST**0 . 5*RSQ(  1 ) ) 

WRITE (6, 1 10 IFF. I  ,ELEF(1) 

1 1 0  FORMAT ( IX ,  1 ***FE 1* 1 . E 1 2 . 5  .  IX .  1 ***ELEF ( 1 . FI 2 . 5 ) 

C 

c _ CALCULATE  THE  SECOND  VALUE  OF  VELOCITY  OF  PARTICLE 

A«-VX(im./ST*DELPS 

B=A*A+4 . * ( VX ( 2 ) *D£LPS/ ST - DELPS*ELEF ( 1 )*COST( 1 ) ) 

CaSQRT(B) 

VPSI (2)*( *A+C)/2 . 

D=l./ST*DELET-VY(1) 

E=D-D+4  ■  '••  ( YY ( 2 )  -DELET/ ST+DELET*EL£F ( 1  )* S I NT(  1 ) ) 

F=SQRT(E1 

VETA(2)*(-D+F)/2. 

DELT*0. 1 
AA*PSI(1) 

BB«ETA(1) 

CALL  MOVEA(AA.BB) 

DO  50  1*1. N 

RSQ(  I+l  )-(PS  1(1+1  J^+ETAC  1  +  1  )**2) 

DIST( I+l )*SQRT(RSQ( I+l ) ) 

COST(I+l)*PSI( I+1)/DIST(I+1) 

SINT(I+1)*ETA(I+1)/DIST(I+1) 

c. . . .CALCULATE  THE  LOCAL  REYNOLD  SO  FROM  THE  RELATION  THAT  REYNOLD  NO 
C  OF  FLOW  FIELD  TIMES  LOCAL  VELOCITY 
VXDIF*(VPSI ( I+l ) »VX( I+l ) )**2 
VYD IF- ( VETA ( I+l ) -VY ( 1+ 1) )**2 
RESOLC I + l ) * R ENOF* ( VXD I F+ V YD I F ) ** 0 . 5 

C0EF( I+l)»l .  +0 .  19  ?*RENOL(  I  +  l  )**0 . 63+0 . 00026’-'RES0L(  1  +  1  )**  1  .  38 

H=COEF (I+l ) 'ST 

KH® SQRT ( ST/  COF.F  ( I  + 1 ) ) 

C 

c. . . .COULOMB  FORCE  NUMBER 
c 

FE®QA*QB*HH/  ( RC-MASSP'-CLENG ) 

C 

c. . . .ELECTRIC  FORCE  TERM 
c 

ELEF (  I  +  l  )*FE*CLENG/  (UINK-"--2--HH-  KS0(  I  +  l ) ) 

C.  .  .  .CALCULATE  THE  DISPLACEMENT  OF  THE  PARTICLE 

PS l  ( 1+2 )®VPSI  ( l  +  l  )*DELT-0.5*H*(VPSIU  +  1  )-VX(  l  +  l )  }«t>ELT**2 
S  +PS I  ( I  + 1 )  -0 . 5 '••ELEF  (I  +  l  V-COST  ( I  + 1  r-'UELT"-"--'2 

ETA(  1+2  )=VETA(  1  +  1  )-'DELT*0 . 5-'H>(VETA(  1  + 1  )-VY(  I  +  l )  i*DELT**2 
$  +ETA(  I  +  l)+0 . 5  ••'ELEF  (I+l ) ->S  I  NT  (I  +  l  )'‘-'DELT**2 

c _ CALCULATE  THE  VELOCITY  OF  PARTICLE 

VPS  I  ( 1+2  )*VPS  I  ( I  + 1 )  -H  -'  ( VPS  I  ( I  + 1 )  -VX  ( I  + 1 )  r-'DKLT 
5  -ELEF (I  +  l  r  COSK  1-1  J'-'DELT 
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VETA ( I +2 )  =VETA  ( I  +  U  -H*(  VETA  ( I  +  U -VY ( I + n  )*DE  LT 
S  +E  LEK 1 1  ♦  1  >  *  S I  NT  il  ♦  1 )  *1  »E  LT 

WRITE(6.24)I,VPSI(I).I.VETA(I) 

24  FORMAT  l  /  IX ,  *  VPS  I  (  '  ,12.')-'  ,E12 .o ,3N, '  VETA(  '  ,  13, '  )«'  ,EL26) 
SQ=SQRT(PSI ( I )**2+ETA( I )**2) 

EPS*ABSfSQW. 

RATIO-RP/RC 

IF(  EPS .  LE  .  RATIO  GO  TO  20 
CA LL  DASHA i  PS  I  <J  ) . ETA UUI 
WRITE(6,22H.PSUn  .  I  ,ETA(H 

22  FORMATl/  IN, '  PSI  ( '  .13,' )='  ,  E 1 2 . 6 , 3N . '  ETA  (  '  .13,')“'  ,f.l2.6) 

GO  TO  50 

20  CALL  ANMODE 

WRITE (6, 2 1)I,PSI(I),I .ETA( I ) 

21  FORMAT  i  /  IN, '  THE  CRITICAL  IT:  ’/lN.TSlC  .13.'  )=’  .E12.6.3X, 

S  ‘ETAi 1 . 13. ' )=' .E12.6) 

GO  TO  55 
50  CONTINUE 
55  RETURN 
END 

C  *•>  ’.'n’rVr  vrYn*;  '.'rCrtrfruirit 
Q  Vf  V.*  -.V  V;  i:  V.*Vr  V;  Vr  '.'"..T;  *  V; 

c _ DRAW  A  CIRCLE  DIVIDED  BY  60  POINTS  WITH  CENTER  0(0.  .0.  ) , PIANETF.R 

c _ CD*  1.  IN  DIMENSIONLESS  WAY 

SUBROUTINE  CIRCLE 
CALL  POINTACO . 0,0.01 
A*2. *3. 1415926/00. 

I1=100*COS(A) 

I2*100*SIN(A) 

CALL  MOVREL(I 1,12) 

DO  5  1=2. ol 
6= A*  I 

IN*100*COS( B) 

I  Y=100,vSINi  B) 

IA=IN-.I1 
I B=IY- 12 
1 1=IX 

::*iy 

CALL  DRWRELdA,  TB> 

5  CONTINUE 
RETURN 
END 
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APPENDIX  -  B 


VISCOUS  FLOW  PAST  AN  INFINITE  CYLINDER 


Schlamp  (1977)  assumed  a  3teady-state  viscous  flow  normal  to 
an  infinitely  long  cylinder  with  the  free  stream  velocity  of 
the  fluid.  The  Stokes-Navier  and  continuity  equations  are  then 


tJ*Vu  -  -VP/p+vV2u  ,  (2-1) 

V  *u  =  0  .  (2-2) 


Considering  well-known  vector  relations  one  can  deduce  from  Eqs. 
(2-1)  and  (2-2) 


Vx(  vV 


2-*- . 
u , 


-Vx  (ux(Vxu)} 


(2-3) 


Velocities  and  other  parameters  have  in  two  dimensional,  polar 
coordinates  the  form 

u*e2xV4i  ;  ur  «  -1/r  (3^/36 )  ;  5  *  7xu  »  ez72 \>  ;  uQ  -  3^/ar  ,  (2-4) 

where  e  is  unit  vector  in  2  direction,  ^  the  stream  function  and 
^  the  vorticity.  We  now  nondimensionalize  the  relationships  in 
the  usual  way 


r'=r/Rc,u«u/uoo,^'*^/u00Rc,NRe*2rpu00/v,c’=Rc/u00  .  (2-5) 

From  Eq.  (2-3)  one  then  finds,  after  dropping  the  primes  for 
convenience , 

V2(7xu)*0.5NRcVX{uX( VXy)  ,  ( 2-S ) 

and  with  equation  (2-4) 

7  (  7xu  )  *0 . 5NRe  7X  { (ezx  V$)x  (  Vx(  ozx  7<p)  )  }  .  (2-7) 

Using  standard  vector  operations,  one  may  write  the 
nond imens ional  Navier  Stokes  equations  as 
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V44/*0.5NRe{  0^/3  r)  o  (7  24/ )/3Q  ] -<3^/39  )  ( (a  (?  >  )/3  r]} /r  .  (2-8) 


More  conveniently,  equation  (2-8)  may  be  split  in  two 
second-order  differential  equations  by  introducing  the  vorticity 
5 .  This  gives 

V2C»0.5NRe{  ( 3 4» /3  r)  (3  ;/39  )- (3^/39  )  (85/3  r)  , 

with  5  =  v2^»  where  7  2  *  3 2/3 r2+l/r 2 ( 3 2/30 2 ) +i( 3 /s r )  .  (2-9) 

r 

Near  the  cylinder  surface  at  small  Re,  a  small  radial  step  size 
is  needed  where  the  stream  function  and  the  vorticity  vary  most 
rapidly.  The  customary  procedure  with  r*ez  was  used  for  a 
constant  step  size  Az.  Equation  (2-9)  then  becomes 

V2;*  0.5Nre{  (3^/32)05/39  )-(3t|;/39  )(3^/3z)  ,  (2-10) 

with  c  ■  e22V2<{>,  where  72  *  32/3z2  +  32/302  . 


Boundary  conditions  are  along  the  axis  of  symmetry:  9*0,  rr,  i(j=0, 
5*0.  On  cylinder  surface:  z*0,  ij/*Q,  C*V2i(i,  3^/3z*0.  On 
cylindrical  boundary  concentric  with  but  remote  from  cylinder 
surface:  2*2,*#  '(’■e2*  sine,  5*0. 

The  flow  pattern  about  an  infinite  length  cylinder  is 
displayed  in  Figure  A. 2-1,  which  shows  that  considerable  fore-aft 
asymmetry  is  already  apparent  for  NRe  =  1.  At  NRf3  *  20,  a 
standing  eddy  of  sizeable  d imens ionsnas  developed  at  the  reverse 
side  of  the  cylinder  (Schlamp,  1977). 


Fig.  B-Ia.  NRe  *  1.0. 
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Fig.  B-Ib.  NRe  «  2.0. 

All  calculations  used  for  the  flow  field  about  an  infinite 
cylinder  were  based  on  the  Ph.D.  thesis  by  Dr.  R.J.  Schlamp 
(UCLA,  1977)  and  on  his  conversion  of  Navier-Stokes  equations  in 
finite  difference  equation.  The  applied  technique  can  be 
summarized  as  follows:  The  first  and  second  derivations  of  basic 
functions  are 


3*/3z  l(IfJ)=U(I,J  +  l)-iJj(I,J-l)}/2A  ,  (2-11) 

V2S/3r2  (I^J)={^(I+1,J)-25(I,J)+;(I-1,J)}/B2  .  (2-12) 


One  can  rewrite  the  equation  (2-9)  using  finite  difference 
formulas  to  obtain 


{  (*(  I,  J+l)-*(  I,  J-l)  )/2A)  (U(I  +  l,J)-UI-l,J))/2B)-(U'(I  +  l,J) 
-*(1-1, J)  J/2B]  (U(I,J  +  1)-  (  I ,  J-l )  ]/2A]  )  = 

2( ( K( I, J+l )-2C ( I , J )+G( I, J-l ) ] /A2] +  [  K( 1+1, J )-2^( I, J) 

+  C(I-1,J))/B2]  }/N'Ro  .  (2-13) 
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This  equation  can  then  be  solved  for  C(I,J)  to  obtain 

C(n  +  1)  (I,J)*[(o.5A2B2/(A2+32)]  U(n)  (I,J+1)  [l/A2+NRe(ti»(n) 

(  I  +  l,J)-4'(n)  (I“l,  J)  1/8A8]  +;(n)  (  I, J-l)  tl/A2-NRe(4»(n)  (I+1,J) 

-4»(n)  (I-L,J)  1/8AB1  +C(n)  (  I  +  1,J)  tl/B2-NRe  [  (4>(n)  (  I,J  +  l)-4,(n) 

(  I,J-1)]8AB]  +  5(n)  (I,  J-l)  [  l/B2-MRe  (  (<|>(n)  (I,J  +  l)-<|/(n) 

( I, J-l) ]/8AB] }  ,  (2-14) 

where  n  and  n+1  introduce  the  concept  of  an  iteration  step. 
Similarly,  if  the  vorticity  equation  (2-9)  is  expressed  using 
finite  difference  formulas  and  solved  for  t|>(i,j),  one  obtains  for 
the  outlined  boundary  conditions 

i|;(n+1)  (I,J)»[0.5A2B2/(A2+B2H  {  (4»(n)(I,J  +  l)-iMI,J-l)l 

/A2]  +  (  [<J>(n)  <I+l,J)+*(n)  (I-l,J)]/B2]-exp[2z(  j  )]  C(n)  (I,J)>  .  (2-15) 

U3ing  equation  (2-10)  one  finds  for  the  side  boundaries 

1  =  1,  M+l  ;  iMl,J)-0,  C(1,J)=0  , 

and  for  the  cylinder  surface 

J-l  ;  iH  1 , 1 )  =0  ,  ;(n+1)(I,l)=0.5A2{84»(n)(I,2)-ti»(n>(I,3)}  . 

For  the  cylindrical  boundary  remote  from  the  collector 
J  =  N;  ii)(I,M)=exp  z(N)sin(I)  ,^(I,:i)=0  . 
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APPENDIX  C 


TWO  HYDRODYNAMIC ALLY  INTERACTING  CYLINDERS 

The  flowfield  between  two  identical  interacting  cylinders 
can  be  described  following  the  scheme  in  Fig.  C-I 
Elementary trigonometric  operations  lead  to  the  deduction  of  the 
velocity  components  in  the  radial  and  tangential  direction  around 
the  reference  cylinder 


Fig.  C-I 

R2  R2 

U  .  j  (i— _)  cose  -  U  — -r  cos  4; 

r  “  x'* 

U  ■  -U  (1+-SU)  Sine  -  U  -A?  sin<J> 
0  ®  r'2 


(  3-1 


where 


r  •  2  s  r  2  +g2  »  2  rs  sin  9.  From  (3-1)  one  obtains 


ur  ■  U  (1- 
r  00 


00  i*  /x2+(8.y)- 


(  3-3) 


and  in  a  similar  way  from  (3-2)  results 


X 


(3-4) 


J? 


u  + 
x 


+y 


/x 


r^uy 

+y 


For  the  numerical  calculation  of  particle  deposition  on 
cylinders  the  flowfield  coordinates  were  converted  into  a 
nondimensional  cartesian  system  with  components 
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1+- 


-2x  y 


(y 


u 


,  ' 2 .  '2, 

(x  +y  ) 


-2x  y 

tT4T7  + 


7^7?  «'*+<|-yV 


9 

rx  s 


(3-^) 


(x  +y  }  /x'  2+y '  2  /xki+(|-  y ' )  ^  lx'2x(|- y'  )2)  (3-6) 


The  dimensionless  parameters  are 


The  calculation  of  particle  trajectory  from  the  flowfield 
and  known  physical  properties  of  the  particle  is  identical  to  the 
procedure  described  in  Appendix  A. 
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APPENDIX  D 


BROWN'S  MODEL  FOR  MANY-FIBER  VISCOUS  FLOW’ 


The  many-fiber  viacou9  flow  model  (Brown,  1934)  Is  based  on 
the  Hehmholt2  theorem  of  the  viscous  flow  pattern — which  adjusts 
to  the  lowest  rate  of  dissipation  of  energy— and  on  the 
formulation  of  a  stream  function  for  a  filter  fiber  array  with 
the  aid  of  the  stream  function  with  spacial  periodicity  (Kirsch 
and  Fuchs,  1967) 

*  =  Hy+fJlE  E  a  .  sin  EJZ  cos  (4-1) 

x  x  n-1  k-0  nk  1 

The  stream  function  was  deduced  for  channel  and  staggered  model 
which  are  depicted  in  Fig.  D-I. 
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formula  for  energy  dissipation 
3u  2  3u  2  3  u 
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f  low 
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3u 


3  z 


3y 


+  <TT  +  T?>  +  <-J7T  ♦  -S^)  >  dxdydz 


(4-2) 


3  z 


3y 


3x 


For  the  two-dimensional  case  with  flow  through  an  array  of  fibers 
extended  in  the  z-direction  Eq.  (4-2)  can  be  rewritten  in  terms 
of  stream  function 

2  2  2  2  2 

9  *  ///{(-^-f  +  2-i)  +  4  (- i)  } dxdy  (4-3) 

3y  3x  3x3v 

2  2  2 
3  >l>  3  d  ip 

— *■  — j  - 

After  calculating  the  terms  of  3y  ,  3x  and  3x3y  from  Eq.  (4-1) 
and  substituting  them  into  Eq.  (4-3)  the  rate  of  energy 
dissipation  will  be 
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CD 


**?  oo  ^  OO  00  2  2 

9  *  n  U  2  it4  et3{  Z  2a  2  +  Z  Z  a  2  (Z-  +  \)}  (4-4) 

“  n-1  no  4  n=l  k=l  nk  2, 2  e 

The  values  of  a  .  that  give  the  minimum  energy  dissipation  and 
satisfy  the  boundary  condition  (stream  function  vanish  at  points 
close  together),  are 


V 

s' 


*>  31  niry.  kTrx.  -y, 

y.+2.  Z  Z  a  ,  — r-i  cos  — -  0  =>  a  .  -  - - - 

1  n=l  k=0  nk  1  e  "k  n*y.  K.x 

Isin — peos — - — 

and  the  condition  of  minimum  energy  dissipation  yield 


(4-5) 


OO  oo 


7  -  -  2  2  t.2  2 

YT —  {  S  2aZ  2_  +  I  Z  a  .  (2-  +  2_)  +• 

3ank  n-1  no  I*  n-1  k»l  k  2, 2  e 


(4-6) 
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00  OO 


niry , 


4  kirx . 

+  Z  y.[y.+2.  Z  Z  a  .  sin— 5-=-  cos  — — -1 }  «  0 
i*l  1  1  n-1  k-0  nK  *  e 


In  Eg.  (4-6)  M  are  the  points  close  to  the  cylinder  surface 
xi'yi#  where  the  stream  function  vanishes  and  y,  are  Lagrangean 
undetermined  multipliers.  Differentiation  with  respect  to  each 
a_i,  gives  an  equation  that  involves  no  other  a_.  .  The  resulting 
equation  for  anR  is 


kTTX  . 


M  nwy 

anR  f(n,k)  +  Z  y  .  2,  sin  -  cos  - 

i-1  1  2  e 


with 


(4-7) 
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Substituting  anR  from  (4-5)  into  (4-7)  yields 
M 

Z  y ■  u • •  =  y •  ,  for  1  <  j  <  M 


i  =  l 


i  17 


with 


(  4-8) 
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UiJ  *  n«l  k=0  £(n'k: 
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n^ry  .  kirx  .  kirx  . 

l  i  i 

— j— cos  -  cos  - 


• ( 4-9) 


Equations  (4-8)  are  solved  by  computer  and  the  results  (  i) 
substituted  into  Eq.  (4-7)  for  determining  a^.  .  Once  ank,fs 
known,  the  stream  function  can  be  calculated  rrom  Eq.  (4-1). 

An  example  of  streamline  pattern  around  a  parallel  array  of 
fibers  for  R  ■  2,  4  -  10,  e  ■  10  is  presented  in  Fig.  D-II 

o-a 

o.» 

«a> 

Fig.  D-II  (Brown,  1984) 

and  for  staggered  model  in  Fig.  D— III  (for  R  *  3,  1  *  10,  e  »  4] 

O.M 

0.8 

O-S 

0.2 
0.02 

Fig.  D-III  (Brown,  1984) 

For  staggered  model  a  stream  function 


-  -  .  "  ?  mry  (2K+l)irx  , 

'P  =  u  v+u  ll  Z  ank  Sin  l  COS  - 2e - 

x  x  odd  n=l  k=0  nx 


as  co 


mry  k-nx 

— COS  - 


+  Ux  1  E  ,  „E„  ank  sin  t  cos  — 

even  n*2  k*0 


is  used  in  the  same  way  as  in  the  previous  model,  except  the 
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odd-n  terms  where  the  following  substitution  is  made 


cos 


knx 


cos 


(2k+l)  irx 


k2)2 


fn2  .  (2k+l) 2, 2 
+  - j-J 

i  (2e) ^ 


When  the  flowfield  in  the  many  fiber  array  is  known  then  the 
particle  deposition  on  each  fiber  and  each  row  (plane)  of 
cylinders  perpendicular  to  the  airflow  direction  can  be 
calculated  according  to  the  relationships  outlined  in  Appendix  A 
(marked  in  Fig.  D-III).  The  main  problem  is  still  the  realistic 
description  of  the  changing  flowfields  inside  the  many-fiber 
scavenger — wh ich  is  related  to  the  pressure  drop  across  the 
filter  (Brown,  1984) — and  the  calculation  of  particle 
trajectories  in  the  many  layer  array.  The  model  seems  to  be  more 
suitable  for  thin  scavengers  with  several  planes  of  fibers  behind 
each  other. 
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APPENDIX  -  E 


PHYSICAL  MODEL  FOR  THE  FALLING  SCAVENGER  ZONE 


The  main  assumptions  taken  as  basis  for  the  numerical  model 
describing  the  evolution  of  a  falling  scavenger  zone  are  as 
follows : 

Scavengers  will  be  dispersed  homogeneously  in  a  zone  located 
at  a  specific  level  at  time  t  »  0  (e.g.,  100  m  or  5000  m) .  They 
start  to  settle,  disperse,  collect  aerosol,  and  entrain  the 
environmental  air.  Scavengers  are  supposed  to  assume  within  one 
second  their  constant  settling  (terminal)  velocity,  which  has 
been  justified  by  our  previous  systematic  experiments. 

The  equations  of  air  motion  and  scavenger  concentration  are 
similar  to  those  suggested  by  Clark  and  List  (1971).  We  also 
used  the  mixing  ratio,  q  (ratio  of  the  scavenger  mass  in  kg  to  1 
kg  of  air)  as  a  measure  of  scavenger  concentration.  In  our 
equations  the  isotropic  turbulence  exchange  coefficient,  K  , 
represents  the  interaction  of  the  induced  air  motion  and  trie 
turbulent  atmosphere.  The  basic  equations  of  the  two  dimensional 
model (Podzimek  and  Smith,  1986)  are: 


3 1  3x  dZ  ax  p 


(£)  +  k.  72u  , 
P  t 


(5-1) 


iz  ♦  u  |2  +  „ 

at  3x  3z 


TZ  (?*■  +  K.  V2w  -  gp 


dZ  p 


(5-2) 


+  22L  *  o 

3  X  3  Z 


(5-3) 


for  air  velocity  components. 


If  *  ^  *  h  ^  -  V>  -  0 


(5-4) 


for  the  scavenger  mixing  ratio,  and 

ff  [j  (u2  +  w2)  +  gqz]  dxdz  =  -  ff  gq  Vg  dxdz  - 

-  Kt  ((!£)2  +  (!r)2  +  (^2  +  (fr)2]  dxd2 

(5-5) 

for  equilibrium  of  energy.  The  last  equation  can  be  rewritten  in 
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terms  of  air  kinetic  energy,  K  *  //  l/2(u^  ♦  w^)dxdz,  Qf 
potential  energy,  gO  *  9  ff  q  dxdz,  go7  ■  g  //  q  a  dxdz  and 
energy  dissipation  (the  last  term  in  Eq.  (5-5)1.  Equations  (5-4) 
and  (5-5)  were  also  used  for  calculation  of  the  convective 
velocity  of  the  center  of  scavenger  mass  zone,  V,  ■  (3z/3t  +  Vs), 
assuming  that  no  scavenger  has  reached  the  ground  (Q  *  const.)? 

Velocity  components,  u,  w  were  calculated  from  the  stream 
function,  , 


(5-6) 


and  are  related  to  the  vorticity  function,  5, 


3 u  _  3w_ 
3  z  3x 


(  5-7) 


by  the  Poisson  equation, 


V2  *  =  5 


(5-8) 


Further,  the  vorticity  transport  equation, 


i£ 

3x 


-  -  H  -  -  H  *  «t  $  *  * 9  & 


(5-9) 


was  derived  by  using  Eqs.  (5-6)  and  (5-7)  to  simplify  Eqs.  (5-1) 
and  (5-2). 

The  initial  and  boundary  conditions  are  apparent  in  Fig.  E-I 
and  were  adapted  to  the  goals  of  specific  numerical  calculations. 
For  each  study  the  size  of  the  velocity  fields  to  be  studied  was 
selected  to  be  (L,  H)  with  a  mesh  spacing  Ax  *  Az  =  A,  and  the 
initial  release  height  of  the  scavenger  zone  to  be  HB  (the  bottom 
of  the  scavenger  zone).  Air  and  scavengers  are  initially  at 


rest.  Then  the  following  boundary 
each  wall  of  the  field  (denoted  by 

conditions 
numbers  in 

were 
Fig . 
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(still  air  boundary):  w 

a  0 ;  u  *  0 , 
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0 ;  (  »  0. 

(5-10) 
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The  following  field  dimensions  (Fig.  E— X)  and  parameters 
were  selected  for  special  investigations  in  the  homogeneous  calm 
atmosphere: 

a)  The  impact  of  scavenger  mixing  ratio,  q,  on  the  induced 
downdraft  velocity,  w,  and  maximum  width  of  the  falling  scavenge 


q  »  0.01  g  sc./g  air  and  0.04 
H  ■  10  km;  Vg 

L  s  10  km; 

HR  a  2  km; 

Zone  dimension:  1  km  x  1  km; 


g  sc./g  air; 

*  1.50  m/s; 

»  1000  m2/s; 

mesh:  20  x  20. 


b)  The  effect  of  changing  the  turbulent 
coefficient,  on  the  downdraft  velocity,  w 
the  scavenger  zone: 


exchange 
and  dispersion 


of 


q  »  0.01  g  sc./g  air; 

H  ■  10  km;  Vg  ■ 

L  =  10  km;  » 

Hr  »  5  km; 

Zone  dimension:  2  km  x 


6.0  m/s; 

1 ; 10; 100; 2 50; 500; 1000; 2500  m2/s. 
2  km;  mesh:  20  x  20. 
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c)  The  influence  of  the  initial  shape  of  the  falling 
scavenger  zone  on  the  induced  velocities,  u,  w,  and  the 
dispersion  of  the  scavengers: 


q  ■  0.01  g  sc./g  air; 

H  -  3  kin;  Vg  -  1.73;1.50;1.21  ra/s ; 

L  -  3  kin;  Kfc  «  1000  mVS ; 

HR  «  1  km  mesh:  30  x  30. 

Zone  dimensions:  400  m  x  400  m,  320  m  x  S00  m, 

500  m  x  320  m,  800  m  x  200  m. 
d)  The  effect  of  the  low  release  height,  H_,  and  of  the 
turbulent  exchange  coefficient,  Kfc,  on  the  evolution  of  the 
downdraft,  w: 


q  «  0.01  g  sc./g  air; 

H  -  300  m;  Vg  «  1.73  m/s; 

L  «  600  m;  «  10  and  1000  m3/3. 

Hr  ■  100  m; 

Zone  dimension:  100  m  x  100  m;  mesh  12  x  24. 

NUMERICAL  MODEL  CRITERIA 

Numerical  solutions  of  the  above  mentioned  equal  Ions  are 
based  on  the  criteria  as  follows: 

Stability  equation  for  determining  (P'ache,  1932.) 

A2  2 

At  <_  — ,  where  a  ■  2216  m  /s. 


Finite  difference  equation  for  the  vorticity  transport  Eq. 
(5-9)  is  using  forward  time,  central  spacing  differencing  in 
fully  explicit  form 


At 
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K 


71  (ei+l,j+^i-l.j+5if jj.+5i.j-l’45i,j>  +  {qi+l,j“qi-l, * 


(5-11) 


The  calculation  of  stream  function,  <J/ ,  is  based  on  the 
successive  over-relaxation  equation  (Frankel,  1950;  Young,  1954) 


„k+l  , k  .  ,  ,k  ,k  .  ,k  »2,  \ 

*i,j  -  n.j +  -  *i,j> 


(5-12) 


with  an  optimal  over-relaxation  factor  approximated  (Frankel, 
1950)  by 

i-yJEK]  .  5  ,  ^coa^)  4-  cos^)^ 


u). 


2  (: 


I,J  are  numbers  of  nodes  in  x-  and  z-direction. 

Scavenger  distribution  and  deformation  of  the  scavenger  20ne 
with  the  time  was  calculated  under  the  assumption  that  regions  of 
scavengers  were  represented  by  a  tracer  particles.  These  were 
expressed  by  the  mixing  ratio,  q  (calculated  for  the  air  mass  in 
a  specific  cell),  divided  by  the  initial  number  of  tracers  in  the 
cell.  At  the  beginning  of  each  time  step,  all  tracers  within  the 
control  zone  of  each  node  are  summed  up  to  calculate  the  new 
mixing  ratio  for  each  node. 

Extending  the  above  mentioned  procedure  into  the  evolution 
of  scavenging  cloud  at  horizontal  wind  velocity  the  following 
boundary  conditions  and  scheme  (in  Fig.  E-II)  were  accepted: 


z- 


I«0 


2 


X.  t.  JO  fcm 


Fig.  E-II 
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At  the  boundary  1 


u  *  U  t  £  ■  0  ;  •  0  , 

w  a  X 

at  the  boundary  2 

w  ■  0  ;  £w  «  ?w+1  (slip  at  the  wall)  ?  «  0  , 

at  the  boundaryy  3 

£  »  0  ;  |>£-  *  0  and 
at  the  boundary  4 

w*0;£»0;i|/a  UQH  (from  u  *  • 
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PROGRAM  ZONE; 

Thu  prograi  ms  d*v*lop*d  for  th*  Graduat*  C*nt *r  for  Cloud  Physics 
R***arcft  at  th*  Univ«r*ity  of  Nissouri-Rolla  und*r  th*  direction  of 
Or.  Jos*f  PooiiMk. 


Scav«ng«r  Particl*  Zon*  r*l»*s*  in  a  Horuontal  Kind 
by  fnoMt  £.  Saith  Jr.  1963  *  1386 

This  progra*  is  a  2-0  finit*  diff*r*nc#  *od*l 

> 

<*R+> 

{*1  GRAPH.P  > 

(II  SMLZXV 
TYPE 

MIX  P  *  "MIX  OAT  ; 

NODE  P  *  "MODE  OAT  ; 

VEL  D  «  "VEL  DfiT  : 

PART  P  3  "PART  OAT  ' 

MIX  SAT  *  AARA?t-l..’&01  OF  ARRAYC-1..201  OF  REAL  ; 

WJOE  OAT  *  ARRAY!*!. .601  OF  ARRAYC-1..201  OF  ARRAYU..21  OF  REhl  : 

VEL  SAT  »  ARRAYt-2. . 1201  OF  AftRAYC-2. . 401  OF  REAL  ; 

PART  OAT  *  ARRAYC1..8001  OF  REAL  ; 

CENTROID  »  ARRAY CO.. SOI  OF  REAL  ; 

VAR 

OUT  ;  TEXT  ; 

MIX.  VIS  :  MIX  P  :  C  MIX  «  lining  ratio  array,  VIS  «  visco*ity  array  > 

VORT,  SAI  ;  NOflE  P  ;  (  VOflT  «  vorticity  array,  SAI  »  strra*  function  array 
IL  U  :  VEL  P  ;  <*U  3  X-dir.  velocity  array,  U  »  Z-dir.  velocity  array  > 

XPOS,  ZPOS-:  PART  P  ;(  XPOS  »  X-coord.  array  of  tr*c*r  particle. 

ZPOS  ■  Z-coord.  array  of  trac*r  particle  } 
ncm,  ZlCAN  :  CENTROID  5  (  BCAN,  ZHEAN  3  coordinate  of  c*ntroiO  of  ton*  > 
MEW,  OLD.  tea,  I.  J,  K,  T  :  INTEGER  ; 

XlTK  Fij«.  INC.  CaKT  :  INTEGER  ;  ^  * 

DELTAT  OELTitt,  VEL,  POS,  ADD,  *IWT  :  REAL  j  (  DELTAT  3  tia.  st*p, 

OELTAX  »  inn  spacing  > 

TIME,  X,  Z  :  REAL  ; 


> 


Initialization  of  array* 


PROCEDURE  INITIVAR  MIX  i  MIX  P;VAR  VORT,S»I  :  NOOE.PjVAR  U.U  :  VEL.P; 
OELTAX  :  REAL!; 

t,  J  :  INTEGER  ; 

Uo  :  REAL  ;  (  Uo  3  antranc*  horizontal  velocity  > 

BEGIN 

Uo  ;«  10.00; 

FOR  J:*0  TO  20  OO 
BEGIN 

FOR  1:*0  TO  60  10 
BEGIN 

NtX"CI, J3  :3  0.0; 

VORT‘Ct.J.11  :3  0.0; 

V0RT"CI.J,21  :3  0.0; 

IF  J»0  THEN 
BEGIN 

SATCI,  J,  11  :=  0.0; 

SAI"::, J, 21  ;s  0.0: 

Offl; 

END; 

END; 

FOR  J:3l  TO  20  OO 
BEGIN 

SAI"CO, J,  1]  :*  SAI"!?,  J*l,  llHIo; 

SAi"CO, J,21  :=  SAI"C0, J, ll; 

FOR  K;*l  TO  2  00 
BEGIN 

FOR  I:*I  TO  60  DO  SAI"C0,J,K3  :=  SAI*C0,J,ll; 

EM); 

END; 
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END; 

FOR  I:=0  TO  60  DO 
BEGIN 

SflIACI, 20,11  :*  SAIAt0,20, 1); 

sni'n,2o,2i  :« sAtAto,20,n; 

EM); 

FOR  l:**2  TO  120  CO 
BEGIN 

FOR  J:*-2  TO  *0  DO 
BEGIN 

Uo: 

0.6; 


U*CI,J) 
•rtt.Jl 
END; 

END; 

END; 

( 


Display  of  Array  Valuas 


) 

PROCEDURE  DISPLAY  I  M  :  NODE  P  ;  t€W  :  INTEGER  ); 

BEGIN 

FOR  J:=20  DOUNTO  0  00 
BEGIN 

GOTOXY  (2,22-J); 

FOR  I: *53  TO  60  DO  WRITE  (  J,NEWU3 j9:3  >; 

EM); 

END; 

< 

Storage  of  Array  Valua* 

) 

PROCEDURE  PRNATRIX  (  N  :  NODE  P  ;  NEW  :  INTEGER  > ; 

VAR 

I.  J,  K  t  INTEGER  ; 

BEGIN 

FOR  K:=0  TO  2  00 
BEGIN 

f  ..  T  |  k< 

FOR  I:*)  l6  GO  00  WRITELN  (  OUT,  rtl, J.NEW3:lO:6,  rCI,K.M€W):10:6, 
r C I ,  Jr2, NEW) :  10:6,  NA C I ,  J*3. NEW) :  10 1 6. 
rd,J*,f€W):10:6,  rCI,J*5, NEW):  10:6. 
NACI,J*i,  NEW):  10:6  !; 

EN0: 

END: 

{ 

Storage  of  Array  Valua* 

\ 

PROCEDURE  PRNATRt:  <  N2  :  VEL  P  ) ; 

VAR 

I,  J,  *  :  INTEGER  ; 

BEGIN 

FGR  K: =0  TC  2  DO 
BEGIN 

J  :s  7  »  K; 

FOR  l:=0  TO  iO  DO  WRITELN  !  CUT,  H2*t2»I,2*J):  10:6.  N2'  :2*I.2«J»2):10:6, 
«2AC2*I,2*a>*);lO:6.  N2'C2»!.i*M3:10:6, 
N2';2*l,2»Jr«) :  10:6,  N2~:2«I.2*.M01 : 10:6, 

M2" '2*1 , 2»J*12) :  10 : 6  ); 

ENO; 

END; 

{ 

Storage  of  Array  Value* 

> 

PROCEDURE  PRNATRU  '  1C  :  NIX  P  ); 

VAR 

I,  J,  K  :  INTEGER  ; 
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8E6IN 

FOR  K:«0  TO  2  CO 
BEGIN 

J  :•  7  *  K: 

FOB  I:«0  TO  60  DO  UBITELN  (  OUT,  H3Atl, J3:10:6.  IC'CI, >11:10:6. 

N3ACI,>23:10:6,  M3MI, >33: 10:6, «ACt, >63: 10:6, 
NTtl,  >51:10:6,  M''tI,J>63:10:6  ): 

END; 

90; 


( 


Storage  of  Particle  Zone  Outline 


PBOCEDUBE  POZOTC  (  XP08,  IPOS  :  PART  P  ); 
VAB 

J,  K  :  INTEGER  ; 

OX,  DZ  :  REAL  ; 

BEGIN 

URJTELN  (  OUT, ’PC  *.?  POSITIONS’  ); 
FOR  Hi *761  TO  V> 

BEGIN 

D1  :*  XP0SA£K3  /  500.0  ; 

OZ  i*  ZPQSACK3  /  500.0  ; 

UBITELN  (  OUT,  01:10:4,  0Z:10:4  ); 
90; 

FOR  J:*19  OOUNTO  2  00 
BEGIN 

K  :•  40  •  J; 

OX  :•  XPOGAOU  /  500.0  ; 

OZ  :*  ZP08*tX3  /  500.0  | 

UBITELN  (  OUT,  0XU0:4,  OZ:iO:4  1; 
90; 

FOR  K:*40  OOUNTO  1  00 
BEGIN 

OX  i*  XPGBAtK3  /  500.0  | 

OZ  :•  ZP08AW3  /  500.0  ; 

UBITELN  i  OUT,  01:10:4,  0Z:10t4  ); 

URiIeLN  (  OUT  ): 

FOB  Ji«2  TO  19  jjO 

BEGIN 

B  :•  (  40  #  1  )  •  39  ; 

OX  s*  XP0SCK3  /  500.0  ; 

OZ  :«  ZPOS'C/.J  /  500.0  ; 

UBITELN  1  OUT,  0X:10:4,  0Z:10:4  1; 
90; 

END; 


Sregfcicai  Duplay  of  Panicle  Zone 

REALiVAR  »fCAN,:«£AN  :  CENTROID; 


PROCEOU*  VIEW1XP0S,  ZPOS  i  PART  PjOELTAX  : 

T  :  INTEGER  );  * 

VAR 

K,  IX,  ZZ  i  INTEGER  ; 

RJLT  i  REAL  ; 

BEGIN 
H1BO  j 

HIRESCOLOR  (9); 

OBAU  (0,0.599.5,1): 

OBAU  <0, 1W.599, 199,  Di 
ORAU  (0,0,0.199.1): 

OBAU  (599,173,573.6.1); 

i*  10. v  /  OEL'PX  ; 
iXCBNCT]  :*  0.0  ; 

ZfCBN'Tl  :•  ■*.  0  j 
FOB  rf.M  '0  800  fiC 
fCGlN 

IF  IP0rACK3  (  0.0  n®  1P08ACKJ  J‘  o.v 
IF  (POS'CKl  )  J0O00.0  THEN  iP0SACK3  I* 


ioooo.o 
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.  v,  /,  f. 


IF  ZPOSAOG  <  0.0  T>€N  ZPO$ACK)  :«  0.0  ; 
XX  :•  ROUND  (  XPOSACK]  *  MULT  )  -  l  ; 

ZZ  :*  199  -  ROUND  (  ZPOS'CKl  •  MULT  ); 
Pun  UX,  ZZ,  1) : 

XlCflf  tTl  :»  XXCflNITJ  ♦  XPOSAtKJ  ; 
ZHEANITT  :»  ZNEANCT]  ♦  ZPOSACIO  ; 

END: 

Oewm  i»  XNEANCT)  /  (  DELTAX  *  800.0  >; 
DCANCT3  :■  ZWAHCT1  /  <  DELTAX  •  800.0  ); 
SOTOXY  (65,23)] 

WRITE  (  XMEANCil:tO:3  ); 

SOTOXY  (65, 24): 

WRITE:  (  D€ANCtlilOt3  ); 


Subroutine  for  interpolation  betueen  node  values 


< 

) 

PROCEDURE  INTERP  XiVAR  VAL  PAW.  PI, VAL  P2,VAL  P3,VAL  PS,I,Z  :  REAL ; 

‘  XI, Z1  :  INTEGER  t; 

VRR 

CONST  A,  VPL  PL,  VAL  PR  :  REAL  ; 

BEGIN 

CONSTA  :■  Z  -  Zl] 

VAL  PL  :«  (  CONSTA  ♦  (  VAL  P2  -  VAL  Pi  i >  ♦  VAL  PI  ; 

VAL’PR  :«  (  CONSTA  ♦  (  VAL"Pt  -  VAL'P2  »)  ♦  VAL‘P3  ; 

VAL'P  :«  U  X  -  XI  )  ♦  (  VfiL  PR  -  VflL  PL  ))  ♦  VAL  PL  ; 

f'  ... 

Subroutine  for  calculation  of  nam  eixing  ratios 

PROCEDURE  NEMIXiVAA  NIX  ;  NIX  P;U,U  ;  VEL  PjVXW  XPOS,2POS  s  PART  P; 

DELTAT,  DELTAX;  VEL  :  REAL") ; 

VAR 

I,  J,  K,  XI,  Zl,  XX.  ZZ,  FLAG  :  INTEGER  ; 

X.  Z.  U  POINT,  U  POINT,  X  TRAVEL  i  REAL  ; 

Z  TRAVEL,  ADO  >  SEAL  ; 

BESlN 

FOR  It^)  TO  60  DO  <  Reset  aetrn  to  0  > 

BEGIN 

FOR  Ji-0  TO  20  DO  NIXACI,J]  !•  0.0} 

EM); 

FOR  K:«l  TO  800  DO 
BEGIN 

X  j«  2.0  •  XPOS  CK3  /  06LTAX  ; 

Z  :•  2.0  •  ZPOS'IKJ  /  OELTAX  ; 

XI  i«  TRUNC  (X); 

Zl  i«  TRlK  (Z){ 

<  Find  velocity  components  at  tracer  locatuns 
INTERP.tlU.POlNT.U-ai.Zll.'J-tXl.ZlMJ.U-CXtM.ZIJ.U'tXlM.ZMj.X.Z.Xl.ZlM 
IMTERP.*<«.POIMT,irCXl,Zl),WA«l,ZmJ,U  IXIM.ZII.ITIXI*:, ZlM), X. Z, XI,  ZS)  i 

<  Determine  distance  of  truer  travel  dun'.;  t;ue  step  > 

X  TRAVEL  i«  U  PCIWT  *  DELTAT  t 

Z, TRAVEL  :•  (’U.POINT  -  'VEL  >  ♦  OIL  TAT  i 
<*  *  Calculate  new  tracer  scsitions  ) 

XPOS*  UU  i •  XPOS* IK)  ♦  X  TRAVEL  ; 

ZP08A[KI  i«  ZPOrtK)  ♦  Z 'TRAVEL  | 

<  Tracers  eust  be  kept  in  ooundary  area 
IF  XPQ6‘IK)  (  g.u  THEN  XPQ6a[K;  i«  0.0  : 

IF  XPOS'IKI  )  WOOO.O  T)€N  XPOS'IKJ  i»  3>000.0  } 

IF  ZPOS'IKI  I  y.O  TiCN  ZPQ6ACX.3  :•  0.0  i 

<  Determination  of  new  uuino  ratios  at  eacr.  n-rte 
U  :•  ROUM)  <(  X FOB* IK]  I  /  DELTAS  )| 

Zl  i»  ROM)  ((  ZP09AM)  )  /  DEL  TAX  >; 

FLAB  i»  0  : 

ADD  i»  O.OuOk  ; 

IF  FRAC  !  XPOSMKJ  /  DELTAX  >  «  0,5  THEN 
BEGIN 

ADO  !•  ADD  *  O.S  ; 

FLAG  i*  FLAG  ♦  l  i 


DC} 

IF  F 


RAC  (  ZPOS'lX]  /  DEL  TAX  >  ■  0,5  T)CN 
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BEGIN 

000  :«  000  *  0.5  ; 

FLOG  : »  FLOS  ♦  2  ; 

EM); 

IF  FLAG  «  0  no  :*  NIXACXI,Z13  *  000  i 

IF  FLOG  «  1  HO 

BC6IN 

NIXACXt,Z13  :■  NIXAH1,Z1J  ♦  000  : 

NIXAtXl-l,Z13  :*  NIX'til -l,  Z13  *  ADO  ; 

E»! 

IF  FLOG  ■  i  n€N 
BEGIN 

MXACX1,ZU  :«  NlXAtXl,ZU  ♦  000  s 
WXAtXl,21-U  :«  NIXACXI,Z1-13  ♦  AOO  i 
OOj 

IF  FU»  »  3  no 
BEGIN 

NIXAH1,Z13  :«  NIXACX1,Z13  ♦  000  ; 

NlXACXt-l,Zt)  :«  NlXAIXl-t,Z13  ♦  000  ; 

NIXACX1,Z1-13  :«  NIX*CX1.Z1-13  ♦  ODD  : 

NIXACXl-l, Zl-13  :«  NIXACX1-1,Z1-I]  ♦  ADD  ; 

DID; 

END; 

EM); 

( 

SuOroutina  to  calculata  non  vorticitin  for  tha  «*»h 

> 

PROCEDURE  (CWV0RT  (NIX,  VIS  :  NIX  PiVOA  VOGT  :  NODE  PjU,U  :  VEL  P; 

DELTAT,  DELTA!  :  FE0L;  NEW,  OLD  i  INTEGER) ;  ' 

VAG 

I.  J  :  INTEGER  s 

GRAVITY,  SUN1,  §UK,  SUU,  SUN,  DELTA  s  REAL  ; 

8ESIN 

GRAVITY  :•  0-009B1  ; 

<  Solution  to  vorticity  transport  aquation  ) 

FOR  Ji*t  TO  19  00 
BEGIN 

FOG  Ij«1  TO  59  00 
BEGIN 

<  Equation  brofcan  up  to  utilix*  e©-proct*«©r  ) 

SUU  !•  0.3»GAAVITY*<NIXACIM,  JJ-NIXAII“l,  J3); 

SUC  i«  v6GTA[IM,j!objtVO^AtlM,J,OLD3MJGTACI,>l,OU)3» 
V0GTAII,M.0L03-*IC; 
sue  i«  vis*n.  ;i*sue/o€Lroxi 

SW3  j»  0.5*rC2»I,6*J3t(VOGTAt:,J-l,OL03-VOGTACI,>t.OU):); 

SUN  :•  0.!*i/'C2»I,i#;3»lVO«TACM,J,OL03*VOGTAC;*l,J,OLO:)| 

SUN  t«  SUU*SWe»$UG*SUN| 

DELTA  <>  DELTAT *SWN / DELTA! i 
V0GTAII,  J,  NEUU3  i*  V0GTAtI, J,OLC1*DELTO; 

END; 

EM); 

EN0| 


Subroutina  to  <J»t train*  n*u  «t raaa  function  .ib»»  uaing 
tna  JOG  aatnod  mth  altamattng  taatp  diracticn 

PROCEDURE  fCUSAI  '  VAR  SOI,  VOGT  i  IflDE  P  |  -OUTOX  >  REAL  ;  VAG  r£UU,  NM  i  INTEGER  ) 
VAR  * 

i.  J,  ccurr.  0LD&,  flag  i  integer  « 

DIF,  sum,  sue,  DELTA,  9 CM  j  re&.  ; 

BEGIN 

UCflA  i*  1.151  i  (  OptiauaP  ovar  raiaiatton  factor  > 

COUNT  i<  0  | 

OLSB  !•  2  •  CUB  | 

FLOG  i*  l  i 

WILE  1  COUNT  1  *00  )  AM  (  FLAG  )  0  )  00 
BEGIN 

DIF  i*G.O  ; 

FLAG  j»  0  | 

(  Potlttva  X*4tr.  taaap  > 

IF  NEW!  ■  1  THEN 
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BEGIN 

FOR  JJ«1  TO  13  00 
BEGIN 

SArC60,J,fCUB)  :«  SAIAC39,  J,f!DB]; 

SAI-'CO, J.NEWBl  :•  SAI'Cl,  J,OLDB); 

FOR  S:*L  TO  59  00 
BEGIN 

sun  :*  VORTAa,JlNEUW]*SQR(OB.TAX); 

SUMS  :>  4.0»SArCI,J, OLOB1  : 

SUC  :«  SArCI*l.J,0L0B]*SAra.J*l,aLDB]*SAIACI-ttJ,ieUB3 

♦sa  r  c : ,  j  - 1 ,  fCWB] -su«i -sue ; 

OELTA  :»  0NEGA*SUN2/4. 0: 

SAI'CI, JflEUS3  :»  SAPCI,  ^OLDBIHELTA; 

EM); 

END; 

END 

<  Negative  X-dir.  weep  > 

Fi  g 
BEGIN 

FOR  J:*t  TO  19  DO 
BEGIN 

SAT  CEO,  J,(€W8J  ;«  SAIAl!9,  J.OLOBI; 

SAI-CO, J-NEW83  s»  SAI'Cl, J, OLDS); 

FOR  I:*53  DOUNTO  1  DO 
BEGIN 

SlMl  :«  VOATACI,J,NEaM]*SQR(DELTAX>; 

SUe  :«  4. 0*SAIACI,J1OLDBli 

sue  »«  SAPCI+l,  J, ^fcy8)+SAIAtI, J»l,0L0BI*SAIAtI-l,  J, OLOBI 
♦SAIACI.J-1,NEU8)-Sun-Sue  ; 

DELTA  i-  OtfBMUe/t.Os 
SArcitJ,fEW0I  :•  SAI'Cl, J,0LD8]*DELTA; 

EM); 

EM); 

EM); 

{  Ch*ck  Fop  convergence  > 

FOR  Ji-l  TO  13  OO 
BEGIN 

FOR  !:•!  TO  53  30 
BEGIN 

IF  FLAG  (  1  HEN 
BEGIN 

DELTA  :•  ABSCSAPCL  J.«N8)-SAIAU, J.OLOB)) ; 

IF  DELTA  )  0.01  THEM  FLAG  i«  1  ; 

END; 

EM); 

END: 

(  DI9PLAT(SA:.:CUB);> 
fCWB  ;•  OLDS; 

OLDB  :•  i-OLDB; 

COUNT  :•  COUTH; 

QOTOXY  (24,24); 

URITELNi’SOR  Nethod  took  COUNT: 3,'  Iterttiont* ) : 

END; 

END; 


( 


SuBreutuo  for  eileuletion  of  X-dir.  velocity  field 


PROCEDURE  NEU  0(8A[  >  NODE  ?;VAR  U  :  VEL  F;tfUB  i  INTEGER (DELTAX  :  REAL); 
vaR  " 

I.  J  i  INTEGER  ; 

BEGIN 

FOR  Ji«0  TO  13  00 

KGIN 

FOR  li-l  TO  33  DO  U*C2»I,*»M1  :»  SAI'll, JH.NEWBl-Wl'IIJ.NEUBJi 
EM); 

FOR  Ji-1  TO  13  00  t/U20,2»J1  :■  1TC;13,2»J); 

FOR  I :  *  l  TO  53  OO 

BEGIN 

iri2*i,40)  »•  .n2*i.3i)  i 

FOR  Jft  TO  13  OO  i«  (U*C2»I,2»M)-y'U»t,2»MJ)/2.0 

*U‘C2*I,2»J-U| 

EM); 

FOR  Ji«l  TO  13  00  (Tt-2,20  :«  O.O-tf'  C2,i*J); 
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FOR  J:«0  TO  *0  OO 
BEGIN 

FOR  I:=0  TO  53  OO  ITKMM.Jl  :*  (lTt2*I*2,  J)HTC2*I,  J3>/2.0HT(2M,  J) 
END; 

END; 


Suoroutina  for  calculation  of  Z-dir.  vtlocity  field 


PROCEDURE  NEW  U(SAI  :  NODE  P;VAR  U  :  VEl  P;f€UB  :  INTEGER;0EITRX  :  REAL!: 
'VfiR  ~ 

I,  J  :  INTEGER  ; 

BEGIN 

FOR  I:«0  TO  59  OO 
BEGIN 

FOR  Ji»t  TO  19  OO  ir(£»I»l,2*J)  i«  SAIAtI,  J,«W8]-SArCM,  J,NEUB3; 
END; 

FOR  J:»l  TO  13  M 
BEGIN 

W^tO,2*J3  :»  iril,2»JJ: 
ira20,2»J]  :•  WAC113,2*J1; 

FOR  t:*l  TO  53  00  irt2*I,2*JJ  :«  <trC2»I*l,2«JJ-WAC2*I-l,2*j])/2.0 

♦4rt2*I*l,2*J1; 

irt-2,2^]  :*  NaC2,2*J1; 

END; 

FOR  1:«0  TO  120  00 
BEGIN 

FOR  J:«0  TO  19  00  WACIl2»Jel3  :•  <lTCIl2»J«-2]-«ACI,2«JI)/2.0+«At;.2*J3 
END; 

END; 

(  Calculation  of  vorticity  at  flip  tall  ) 

PROCEDURE  VORTBCIVAR  VORT  :  NODE  P;«WU  :  INTEGER) ; 

VflR  * 

I  :  INTEGER  ; 

BEGIN 

FOR  Is*»  TO  53  00  V0RTACI,0.)«WW  :•  VQKTAU,  l,l«Ml 
END: 

PROCEDURE  SIGNAL  ; 

BEGIN 

SOM  (440); 

DELAY  (300); 

NOGOUNO  ; 

DELAY  (30); 

SOUND  (SOOI; 

OELAY  (500); 

NOSOUND  ; 

DELAY  (30); 

SOUND  (320) ; 

OELAY  (300); 

NQ80UND  ; 

END: 

BEGIN 
)CU  (Nil); 

MEU  (VIS); 

NEM  (VORT); 

«U  <SA1); 

Nft  (U); 

)CU  (U)j 
NEW  (IPOS); 

*CW  (2P0S) ; 

(  Initial iz«ttor>  of  viKcaity  field  ) 

FOR  lt*>  TO  60  OO 
BC9IN 

FOR  Ji«0  *0  10  DO  VIS'CI.JI  )•  1000.0  ( 

00; 

!  Sotting  of  paraMtm  > 

'.EL  i»  6.  .'  ; 

OCLTfiT  i ■  28.0  i 
OELTM  i  •  500.0  i 
OLD  i»  »  i 
NCWN  !•  1  t 
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».EUB  :>  •  ; 

INC  :M  ; 

COUNT  :*  1  ; 

'NIT  (  NIX,  VORT  SAI,  U,  U,  DELTA!  ); 

{  Creation  of  tor*  tracer  particles  > 

FOR  J: -0  TO  13  DO 
BEGIN 

FOR  I:*I  TO  *0  00 
BEGIN 

POS  :«  (  I  ♦  100.0  1  4  9950.0  ; 

K  :*  (  SO  *  J  !  M  ; 

XPQS'CKl  :*  PCS  ; 

END; 

END; 

FOR  l:«0  TO  23  DO 
BEGIN 

FOR  J:*'.  TO  20  00 
BEGIN 

AOS  :»  (  J  ♦  100.0  )  ♦  6950.0  ; 

*  ;«  l  V)  ♦  <  J  -  l  ))  MM; 

ZIOS'OO  :*  POS  ; 

END; 

END; 

{  Calculation  of  initial  suing  ratios  > 

FOR  K:«l  TO  800  00 

BEGIN 

1  :*  XPOS'CK]  /  DELTA!  ; 

’  :•  !P0S*CK1  /  OELTAX  ; 

U  :«  TRUNC  IX); 
il  :«  TRUE  1 2) ; 

FLA6  :*  0  : 

ADD  :*  O.OOO*  : 

IF  FRAC  (X)  *  0.3  THEN 

begin 

ADD  ;*  ADO  ♦  0.3  ; 

FLAG  :  *  rLA6  *•  I  ; 

END; 

IF  FRAC  >.l)  *  0.3  TIEN 
3E3IN 

ADO  :•  ADO  *  0.3  ; 

FLAG  ;•  FLAG  *&  t 
210; 

:?  FLAG  •  0  THEN  nix  '1X1,111  j»  mx  CXI,  112  ♦  ADD  ; 
iF  :LA5  *  1  THEN 

begin 

«ix  :xi, ;n  :» in  cxi.zi:  ♦  ado  ; 

;«  nix'hi-;,:u  ♦  Add  ; 

2*D; 

:?  flag  « i  hen 

BEGIN 

Kirrxuij  :•  wr»i.zu  ♦  add 

«XADXi,2l-tJ  i«  N[X*a;,21-U  »  ADD  ; 

END; 

:f  flag  m  tien 
begin 

iix':xi,:i:  **  iixaixi.:i;  *  add  ; 
in  txK'ij  i»  «x\'xi-i.zn  »  Add  ; 
nx  i<'.,:!-u  ;•  iixAtxi.:i-i:  ♦  add  ; 

-:x  >«  nx'  :x;-i,zi-u  *  Add  i 

END; 

,IE A  :  l?03.  EPOS,  OilTAX,  XfEAN,  ZCAN,  0  >i 
:teppi/.c  loop  through  tiM  / 

FOR  m«1  f0  26  DO 
BEGIN 

T Iff  •  *  OELTAT  *  k  ; 

OLD  )•  J  -  OLD  | 

NEW  ;•  3  -  NEW  j 

:.El*:x  ■  m,  l,  a.  (POfl,  ZW8,  deitat,  delta*,  VEL 
view  !  1P09,  IPOS,  OELTAX,  XKAN,  MAN,  k  ■; 

30T0XV  - TC, 1) ; 
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30TOXY  (76,1); 
kRITE  i  TIME : 4 : l  ); 

NEW0RT  <  Nil,  Via,  VORT,  U,  W,  DELTAT,  CELT  AX,  NEW,  OLD 
VORTBC  (  VORT,  NEVfal  ): 

«USAI  (  SAI,  VORT,  OELTAX,  NEUU,  CUB  > ; 

NEW  U  (  SAI,  U,  NEWS,  DEL TAX  ) ; 

*W/  (  SAI,  U,  fCW,  OELTAX  >; 

{  *  Storing  of  output  > 

IF  INC  *  12  TiCM 
BEGIN 

IF  COUNT  «  1  T>€N  ASSIGN  (  OUT, * 300- ICO'  >; 

IF  COUNT  a  2  THEN  ASSIGN  (  OUT, '600. 100'  ); 

IF  COUNT  *  3  THEN  ASSIGN  :  OUT,  '900.100'  >; 

REWRITE  i  CUT  )j 

WRITELN  (  XT,'  Matrix  ElMtnts  at  \TII€:5:1,’  Socor.as' 
WRITELN  (XT,'  '  ); 

WRITELN  (  XT,'  MIXING  RATIOS’  >; 

PRMATRX3  (  MIX  >; 

WRITELN  i  XT.’  VORTICITIES’  ); 

PfWATRJX  I  VORT.NEWW  /; 

WRITELN  (  XT STREAM  FUNCTIONS'  ); 

PRMATR1X  (  SAI.NEUB  > ; 

WRITELN  (  OUT,1  X-DIRECTION  VELOCITIES’  ); 

PRMATRX2  <  U  ): 

WRITELN  (  OUT/  Z -DIRECTION  VELOCITIES’  >; 

PRMATRX2  (  W  ) ; 

PRItIC  (  XPOS,  ZPOS  ) ; 

INC  S3  0  : 

COUNT  !•  COUNT  ♦  1  | 

CLOSE  <  OUT  ) j 
EM); 

INC  *■  INC  ♦  1  ; 

END; 

ASSIGN  i  OUT,' CENTIF.  .00*  ); 

'-write  (  xt  »: 

-’ITT-',  i  OUT  •  iOfC  CENTROID  (  x, 2  )'  >: 

FTP.  I:")  TO  40  CO  WRITELN  <  OUT,  XMEAMCIS :  10:2, ’  V^ANdl 
CLOSE  XT  i ; 

SIGNk.  ; 

Close  ^  eut  ; 

•OCTOX‘;  .78,24); 

END. 
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